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1.  Statement  of  the  problem  studied 


Objective 

The  objective  of  this  ARO  project  is  to  lay  the  control-theoretic  foundation  of  controlling  a 
hyper-redundant  robotic  arm  or  flexible  robotic  actuators,  often  represented  by  a  distributed 
parameter  system  or  a  partial  differential  equation.  We  establish  a  model-guided  approach  for 
fabrication  and  control  of  flexible  3-D  cellular  actuator  structures  based  on  Electro-Active 
Polymers  (EAPs)  or  other  flexible  materials.  A  model-guided  approach  for  control  of  flexible 
structures  is  epitomized  by  the  Pi’s  PDE  boundary  control  framework  applied  to  articulated 
structures.  The  combined  theoretical  and  experimental  modeling  framework  would  enable  the 
exploitation  and  optimization  of  different  actuator  designs  to  achieve  desirable  dynamic  and 
material  characteristics.  As  far  as  material  characterization  is  concerned,  we  focused  on:  (i) 
multi-physics  modeling  of  EAP  under  combined  electric  and  mechanical  loading  and  (ii)  the  role 
of  geometric  reinforcement  and  stiffeners  on  macroscopic  response  and  force-stroke 
characteristics  of  an  EAP -based  actuator. 


Approach 


•  Derive  a  control-theoretic  framework  of  controlling  a  distributed  parameter  system  (DPS)  or 
partial  differential  equation  (PDE)  using  the  newly  proposed  PDE  boundary  control 
framework  of  Dyadic  Perturbation  Observer  (DPO)  developed  by  the  PI  and  his  research 
team. 

•  Apply  DPO  and  compare  with  the  existing  PDE  control  approaches 

•  Develop  a  simple,  effective  experimental  validation  of  the  otherwise  highly  control  theoretic 
approach. 

•  Devise  a  numerical  framework  that  enables  the  optimization  and  characterization  of  different 
actuator  configurations 

•  Verify  the  material  model  prediction  through  different  combination  of  loading. 

Examine  the  role  of  geometric  constraints  on  the  macroscopic  response  of  the  proposed 
actuators,  to  control  the  force-stroke  characteristics 

•  Devise  a  simplified  analytical  framework  to  describe  the  actuator  response  that  can  be  used 
with  the  control  algorithms. 

Relevance  to  Army 

Compliant  robotic  actuators  provide  the  ability  to  interact  with  dynamically  varying  external 
conditions.  The  realization  of  the  proposed  actuation  and  sensory  framework  will  assist  in 
exploration  of  forbidden  terrains,  monitoring  of  critical  assets,  as  well  as  mitigation  of 
environmental  hazards.  The  Army’s  future  combat  robots  equipped  with  flexible  continuum 
actuators  can  respond  to  combat  needs  in  hostile  environments,  especially  when  integrated  with 
intelligence,  surveillance,  reconnaissance  (ISR)  sensors.  Moreover,  the  proposed  work  can  be 
further  extended  to  computer-controlled  minimally-invasive  surgical  tools  for  the  Army’s  field 
hospitals.  Additionally  actively-controlled  planar  actuators  are  identified  as  enabling 
technologies  for  the  broader  problem  of  developing  a  flexible-winged  micro  aerial  vehicle 
(MAV)  capable  of  agile  flight  in  constrained  environments.  At  the  commercial  and  civilian 
application  level,  we  envision  a  range  of  system-level  applications,  including  (a)  a  generalized 


3D  morphing  structure  for  robotic  grippers,  (b)  flexible  robot  manipulators  without  articulated 
joints  for  greater  maneuverability  in  highly  constrained  spaces  for  endoscopic  surgical 
procedures,  (c)  artificial  muscles  for  biomimetic  locomotion,  swimming,  and  flying,  and  (d) 
next-generation  robotic  manipulators  and  tethers  for  autonomous  docking  in  space  and  tethered 
spacecraft  operations. 

2.  Summary  of  the  most  important  results 

Detailed,  bulleted  accomplishments  for  the  period  should  be  traceable  back  to  approach  and 
objectives  please  be  specific  and  quantitative  where  possible. 

•  We  have  obtained  important  results  on  how  to  synthesize,  analyze,  and  design  control  of 
complex  flexible  structures  by  PDE  boundary  control  formulations  (IEEE  Transactions  on 
Robotics  full-length  paper  published  in  2013).  This  result  is  particularly  useful  when  the  PDE 
model  (combined  hyperbolic  PDEs  of  torsion  and  fourth-order  bending  PDEs)  is  cross- 
coupled  with  the  rigid  body  motion  of  a  platform  represented  by  a  highly  uncertain  nonlinear 
ODE  model. 

•  We  also  derived  and  further  developed  the  newly  proposed  PDE  control  framework  of 
Dyadic  Perturbation  Observer  (DPO)  and  submitted  a  manuscript  to  the  IEEE  Transactions 
on  Automatic  Control. 

•  At  the  same  time,  we  are  making  good  progress  in  applying  a  Central  Pattern  Generator- 
based  control  approach  to  control  and  motion  planning  of  a  flexible  robotic  actuator  that 
might  not  engage  in  oscillatory  motions.  Hence,  this  work  is  different  from  Pi’s  prior  work 
that  has  successfully  applied  CPG-based  control  approaches  to  oscillatory  motions  like 
flapping. 

•  We  developed  mathematical  framework  for  the  electro-mechanical  behavior  of  EAP, 
combining  the  Maxwell-stress  arising  from  applied  electric  field  with  the  non-linear  hyper¬ 
elastic  material  response 

•  We  implemented  the  mathematical  framework  into  a  user-material  subroutine  for  ABAQUS 
commercial  finite  element  software. 

•  We  performed  multi-axial  tensile  testing  with  different  biaxial  loading  ratio  to  calibrate  the 
material  model 

•  We  performed  controlled  actuation  of  pre-stretched  actuators  to  calibrate  the  electro¬ 
mechanical  coupling  and  the  FEM  model  prediction. 

•  We  designed  and  optimized  different  actuator  configuration  capable  of  multiple-force  stroke 
characteristics 

•  We  utilized  segmented  electrodes  to  build  free  standing  planar  actuators  capable  of  3D 
deformation  with  prescribed  force-stroke  characteristics. 

•  We  developed  performance  map  for  different  geometric  stiffeners  layout  to  control  the  force- 
stroke  characteristics  at  constant  overall  actuator  stiffness  and  applied  electric  field. 

•  A  tenfold  increase  in  the  force-stroke  characteristics  is  achieved  from  changing  the  stiffeners 
pitch  and  periodicity. 


•  We  developed  and  verified  simplified  anal5dical  models  based  on  Timoshenko  biomaterial- 
thermostat  theory.  An  experimentally  calibrated  geometric  factor  is  used  to  correct  for  local 
inhomogeneous  deformation. 

Collaborations  and  Technology  Transfer 

•  Invited  talk  at  the  Army  Research  Lab  (Host:  Dr.  Sam  Stanton  &  Harris  Edge)  on 
11/1/2012. 

•  Discussion  with  the  ARL  researchers  on  the  Army  MAST  program. 

•  Data  sharing  with  Pratt  &  Whitney  for  possible  implementation  of  long  slender  internal 
or  external  gripper 

•  Data  sharing  with  Boeing  St.  Lewis  for  possible  utilization  of  soft  actuators 

3.  Bibliography 

Resulting  Journal  Publications  (attached  in  this  report) 

■  A.  A.  Paranjape,  J.  Guan,  S.-J.  Chung,  and  M.  Krstic,  “PDE  Boundary  Control  for  Elexible 
Articulated  Wings  on  a  Robotic  Aircraft,”  IEEE  Transactions  on  Robotics,  vol.  29,  no.  3, 
June  2013,  pp.  625-640. 

■  A.  A.  Paranjape  and  S.-J.  Chung,  “Dyadic  Perturbation  Observer  Framework  for  Control  of  a 
Class  of  Nonlinear  ODE/PDE  Systems,”  IEEE  Transactions  on  Automatic  Control,  under 
review  (submitted  8/2013),  2013. 

■  Lai,  W.;  Bastawros,  A.;  Hong,  W.,  Chung,  S.  2013,  Planar  dielectric  elastomer  actuator 
capable  of  out-of-plane  deformation  without  prestretch,  (In  preparation) 

■  Lai,  W.;  Bastawros,  A.;  Hong,  W.,  Chung,  S.  2013,  Control  of  force-stroke  characteristics  of 
planar  dielectric  elastomer  actuator  (In  preparation) 

Resulting  Conference  Publications 

•  Paranjape,  J.  Guan,  S.-J.  Chung,  M.  Krstic,  “PDE  Boundary  Control  for  Euler-Bernoulli 
Beam  Using  a  Two  Stage  Perturbation  Observer,”  IEEE  Conference  on  Decision  and  Control 
(CDC),  Maui,  HI,  December  2012,  pp.  4442-4448. 

•  A.  Paranjape,  S.-J.  Chung,  and  M.  Krstic,  “PDE  Boundary  Control  for  Flexible  Articulated 
Aircraft  Wings,”  AIAA  Guidance,  Navigation,  and  Control  Conference,  Portland  OR,  August 
2011,  AIAA  201 1-6486. 

•  Lai,  W.;  Bastawros,  A.;  Hong,  W.,  “Out-of-Plane  Motion  of  a  Planar  Dielectric  Elastomer 
Actuator  with  Distributed  Stiffeners,”  Conference  on  Electroactive  Polymer  Actuators  and 
Devices  (EAPAD),  San  Diego,  CA,  MAR  12-15,  2012;  Proceedings  ofSPIE,  Vol.  8340, 
Article  Number:  834011. 

•  W.  Lai,  A.  Bastawros,  W.  Hong,  and  S.-J.  Chung,  “Fabrication  and  Analysis  of  Planar 
Dielectric  Elastomer  Actuators  Capable  of  Complex  3-D  Deformation,”  IEEE  International 
Conference  on  Robotics  and  Automation  (ICRA),  St  Paul,  MN  May  14-18,  2012. 


Graduate  Students  and  Thesis  Titles 

•  Jin}^!  Guan  (M.S.,  graduated  in  May  2012),  Aerospace  Engineering,  UIUC 

M.S.  Thesis:  Design  and  Control  Strategy  of  a  Flexible,  Hyper-Redundant  Robotic  Arm 
using  Electroactive  Dielectric  Polymers 

•  Nitish  Sanghi  (M.S.,  graduated  in  May  2013),  Mechanical  Engineering,  UIUC 
M.S.  Thesis:  Design  of  Flexible  Actuators  using  Electro-Active  Polymers  And  CPG- 
based  Control  Strategies 

•  Michael  Dorothy  (Ph.D.,  Expected  December  2014)  Aerospace  Engineering,  UIUC 
Ph.D.  thesis  in  progress:  Nonlinear  Stability  of  Coupled  Oscillators 

•  Dr.  Aditya  Paranjape  (PhD.,  graduated  in  12/201 1  and  currently  postdoc  at  UIUC) 

Ph.D.  Thesis:  Dynamics  and  Control  of  Robotic  Aircraft  With  Articulated  Wings 

Dr.  Aditya  Paranjape,  who  was  partially  supported  by  this  ARO  grant  will  become  a 
tenure-track  assistant  professor  of  Mechanical  Engineering  at  McGill  University, 
Montreal,  Canada,  starting  in  2014. 

•  William  Lai  (M.S.,  graduated  in  May  2011),  Aerospace  Engineering  (Major:  Engineering 
Mechanics),  Iowa  State  University. 

M.S.  Thesis:  Characteristics  of  dielectric  elastomers  and  fabrication  of  dielectric 
elastomer  actuators  for  artificial  muscle  applications. 

•  William  Lai  (PhD,  Expected  June  2014),  Analysis  and  Eabrications  of  Planar  Dielectric 
Elastomer  Actuators  Capable  of  Complex  3D  Deformation,  Aerospace  Engineering,  Iowa 
State  University 


4.  Appendixes 


IEEE  TRANSACTIONS  ON  ROBOTICS,  VOL.  29,  NO.  3,  JUNE  2013 


625 


PDE  Boundary  Control  for  Flexible  Articulated 
Wings  on  a  Robotic  Aircraft 

Aditya  A.  Paranjape,  Member,  IEEE,  Jinyu  Guan,  Soon-Jo  Chung,  Senior  Member,  IEEE, 

and  Miroslav  Krstic,  Eellow,  IEEE 


Abstract — This  paper  presents  a  boundary  control  formulation 
for  distributed  parameter  systems  described  by  partial  differential 
equations  (PDEs)  and  whose  output  is  given  by  a  spatial  integral 
of  weighted  functions  of  the  state.  This  formulation  is  directly  ap¬ 
plicable  to  the  control  of  small  robotic  aircraft  with  articulated 
flexible  wings,  where  the  output  of  interest  is  the  net  aerodynamic 
force  or  moment.  The  deformation  of  flexible  wings  can  be  con¬ 
trolled  by  actuators  that  are  located  at  the  root  or  the  tip  of  the 
wing.  The  problem  of  designing  a  tracking  controller  for  wing  twist 
is  addressed  using  a  combination  of  PDE  backstepping  for  feed¬ 
back  stabilization  and  feed-forward  trajectory  planning.  We  also 
design  an  adaptive  tracking  controller  for  wing  tip  actuators.  For 
wing  bending,  we  present  a  novel  control  scheme  that  is  based 
on  a  two-stage  perturbation  observer.  A  trajectory  planning-based 
feed-forward  tracker  is  designed  using  only  one  component  of  the 
observer  whose  dynamics  are  homogeneous  and  amenable  to  tra¬ 
jectory  planning.  The  two  components,  put  together,  estimate  the 
external  forces  and  unmodeled  system  dynamics.  The  effectiveness 
of  the  proposed  controllers  for  twist  and  bending  is  demonstrated 
by  simulations.  This  paper  also  reports  experimental  validation  of 
the  perturbation-observer-based  controller  for  beam  bending. 

Index  Terms — Distributed  parameter  systems,  nonlinear  control 
systems,  robot  control,  robot  motion,  unmanned  aerial  vehicles. 

I.  Introduction 

There  is  a  considerable  interest  in  developing  robotic  air¬ 
craft,  which  is  inspired  by  birds  and  bats  [10],  [16],  [20], 
[43]  and  insects  [13]-[15],  [60].  While  insect  wings  can  be  mod¬ 
eled  as  simple  rigid  wings,  both  wing  flexibility  and  wing  artic¬ 
ulation  are  believed  to  play  a  key  role  in  flight  performance  and 
agility  for  bird  and  bat  flight  [43],  [54].  A  new  concept  to  con¬ 
trol  microaerial  vehicles  (MAVs),  which  uses  wing  articulation 
and  is  inspired  by  bird  and  bat  flight,  was  introduced  by  Paran¬ 
jape  and  coauthors  [39],  [40].  The  concept  lends  itself  readily 
to  aeroelastic  tailoring,  which  is  seen  as  an  asset  in  the  develop¬ 
ment  of  agile  MAVs  [26],  [39],  [50],  [55].  Wing  flexibility  not 
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Eig.  1 .  Robotic  bat  testbed,  where  the  control  laws  proposed  in  this  paper  can 
be  tested  [10],  [16]. 

only  improves  aircraft  performance  and  stability  passively,  but 
can  also  be  actuated  actively  for  control  [1],  [39],  [55].  Eig.  1 
shows  a  robotic  bat  testbed,  which  was  developed  at  the  Uni¬ 
versity  of  Illinois  at  Urbana-Champaign  (UIUC)  by  the  authors 
of  [10]  and  [16].  Each  wing  is  actuated  at  the  root,  and  during 
experiments,  the  motion  of  each  wing  as  well  as  the  body  can  be 
measured  in  real  time  for  feedback  control.  It  provides  a  suitable 
platform  to  test  the  controllers  that  are  developed  in  this  paper. 

Broadly  speaking,  continuum  distributed  parameter  systems 
(DPS)  such  as  flexible  wings  can  be  controlled  in  two  ways: 
distributed  control  and  boundary  control.  The  former  approach 
relies  on  a  series  or  a  network  of  actuators  and  sensors  that  are 
distributed  over  the  system.  The  latter  approach  relies  on  actua¬ 
tors  that  are  deployed  only  on  the  boundary  of  the  system  under 
consideration.  Eor  practical  engineering  applications,  the  most 
important  benefit  of  boundary  control  is  the  reduction  in  the 
number  of  actuators.  Despite  a  number  of  bio-inspired  exam¬ 
ples  of  distributed  actuation,  it  is  evident  in  robotic  applications 
that  implementing  a  distributed  actuation  scheme  would  incur 
substantial  penalties  in  weight  as  well  as  costs.  A  distributed  ar¬ 
ray  of  sensors  is  considerably  cheaper  and  less  cumbersome  to 
implement  than  an  array  of  distributed  actuators.  The  benefits  of 
distributed  actuation  can  be  made  up,  in  part,  by  a  combination 
of  good  mechanical  design  and  effective  boundary  control. 

A.  Literature  Review 

Structural  stability  of  the  wing  is  governed  by  its  geome¬ 
try  and  speed.  Moreover,  a  typical  flexible  wing  has  at  least 
two  elastic  degrees  of  freedom  for  deformation — bending  and 
twisting — in  addition  to  the  rigid  rotations  at  the  root.  The  struc¬ 
tural  dynamics  of  the  wing  are  coupled  to  the  rigid  body  dynam¬ 
ics  of  the  aircraft  [39].  It  has  been  demonstrated  in  the  literature 
that  aeroelastic  instabilities  such  as  wing  divergence  and  flutter 
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can  be  mitigated  using  flap-based  effectors  [3],  [46]  or  passive 
energy  sinks  for  flutter  [25]. 

There  is  a  substantial  amount  of  literature  on  boundary  control 
theory  of  partial  differential  equations  (PDEs)  (see  [8],  [9],  [12], 
[30]-[32] ,  [47]  for  material  pertinent  to  this  paper  and  references 
cited  therein).  There  are  two  sets  of  methods  for  boundary  con¬ 
trol  of  PDEs.  The  first  set  consists  of  methods  that  convert  the 
PDEs  into  ordinary  differential  equations  using  approximation 
methods  such  as  those  of  Galerkin  or  Rayleigh-Ritz  [9],  [22], 
or  using  operator  theoretic  tools  [4],  [5],  [11],  [27],  [33],  [34]. 
The  second  set  of  methods  leaves  the  PDEs  intact,  and  uses  a 
“model-following”  approach  as  described  in  a  recent  book  by 
Krstic  and  Smyshlyaev  [30]. 

A  finite- state  approximation  to  a  PDE  may  wrongly  render 
fundamental  system  theoretic  properties  like  controllability  and 
observability  into  functions  of  the  approximation  [8].  Stabil¬ 
ity  analysis  that  is  based  on  a  finite-state  approximation  is 
vulnerable  to  spillover  instabilities  which  arise  due  to  inade¬ 
quate  accounting  of  the  residual  modes  [2],  [37].  On  the  other 
hand,  keeping  PDEs  intact  makes  the  control  law  design  more 
intuitive.  A  PDE-based  approach  has  been  used  for  maneu¬ 
vering  robotic  arms  [30],  [44],  controlling  the  Navier-Stokes 
model  [58],  and  suppressing  vibrations  in  a  flexible  beam  [21]. 
A  gain- scheduling-based  approach  for  nonlinear  PDEs  has  been 
presented  in  [52],  while  Krstic  and  Smyshlyaev  [31]  derived  an 
adaptive  controller  for  parabolic  PDEs. 

There  are  several  PDE  and  ODE  controller  designs  in  the 
literature  to  stabilize  and  control  beam  bending  [17],  [21],  [28], 
[33],  [34],  [52],  [56],  although  the  literature  mostly  concerns 
bending  that  is  encountered  in  loaded  cables,  structural  beams, 
or  robotic  arms.  Elexible  wings  of  robotic  aircraft  are  unique 
because  they  experience  a  closed-loop  interaction  with  the  flow 
field  around  them.  The  loop  closure  happens  due  to  the  flow  field 
itself  being  dependent  on  the  wing  deformation,  which  creates 
a  highly  nonlinear  system  whose  modeling  is  susceptible  to 
significant  structural  and  parametric  uncertainties.  This  neces¬ 
sitates  a  control  design  approach,  where  the  control  parameters 
can  be  tuned  with  very  limited  information  about  the  system, 
and  traditionally  successful  and  applied  methods  for  DPS,  such 
as  [5],  [33],  [34],  need  to  be  replaced  by  a  more  “adaptive” 
approach,  which  can  tune  the  control  parameters  to  appropriate 
functions  of  the  system  parameters. 

B.  Main  Contributions 

This  paper  presents  a  boundary  control  problem  for  wing  twist 
(second-order  hyperbolic  PDE)  as  well  as  bending  (fourth-order 
Euler-Bernoulli  beam  PDE),  which  could  be  applied  to  a  wider 
class  of  PDEs,  including  flexible  aircraft  wings  and  robotic 
arms.  One  highlight  of  our  approach  is  that  it  requires  minimal 
information  about  the  system  dynamics,  and  the  control  laws  are 
themselves  in  quite  a  tractable  form,  which  makes  them  easy  to 
apply  to  a  real  physical  system.  The  stability  of  the  controlled 
system  is  analyzed  rigorously. 

The  main  contributions  of  this  paper  are  as  follows. 

1)  The  PDE  control  objective  that  is  considered  in  this  paper 
is  unique  since  the  controller  minimizes  a  tracking  error  of 


a  time- varying  net  aerodynamic  force  or  moment,  rather 
than  following  a  deformed  wing  shape.  An  abstraction 
of  this  control  objective  is  to  ensure  that  some  weighted 
integral  of  the  twisting  angle  equals  the  prescribed  value. 
A  similar  idea  for  a  more  general  class  of  PDEs  has 
been  explored  in  [8]  and  [9]  in  the  context  of  spatially 
distributed  control. 

2)  We  design  a  tracking  controller  for  wing  twist  based  on 
trajectory-planning  and  a  PDE  backstepping-based  sta¬ 
bilizer.  The  procedure  is  called  backstepping  because  it 
involves  a  Volterra  operator  with  a  lower  triangular  struc¬ 
ture  similar  to  backstepping  transformation  for  ordinary 
differential  equations  [29],  [30],  and  allows  the  controller 
acting  at  the  boundary  to  compensate  for  the  undesired 
(unmatched)  dynamics. 

3)  With  the  integral  of  wing  twist  as  the  output,  and  shear 
force  at  the  wing  tip  as  the  input,  we  show  that  the  twist 
dynamics  have  a  flnite  relative  degree,  thereby  permitting 
a  traditional  ODE-based  approach  to  control  design.  The 
tracking  controller  increases  the  stability  margin  of  the 
twist  dynamics  by  a  factor  of  16.  An  adaptive  controller 
for  a  limited  class  of  parametric  uncertainties  is  also  de¬ 
rived  for  tip-based  actuators  in  general,  and  root-based 
actuators,  when  the  output  is  the  rolling  moment. 

4)  We  present  a  novel  perturbation-observer-based  controller 
to  facilitate  a  trajectory -planning-based  tracking  con¬ 
troller  for  bending.  The  output  of  interest  is  the  displace¬ 
ment  of  the  wing  tip.  The  perturbation  observer,  which 
is  designed  here  uses  adaptation  to  estimate  the  external 
forces  acting  on  the  system.  The  observer  is  split  into  a 
“particular”  component  and  a  “homogeneous”  component 
(the  notions  are  made  more  precise  later).  Since  the  ho¬ 
mogeneous  component  is  stable  and  not  driven  directly 
by  external  feedback,  it  is  simpler  to  design  a  control  law 
for  it.  The  same  control  signal  is  sent  to  the  actual  sys¬ 
tem,  whose  states  then  converge  exponentially  to  a  region 
around  the  observer  states. 

5)  We  experimentally  validate  the  perturbation-observer- 
based  controller  for  a  beam,  which  is  controlled  at  the 
wing  boundary. 

This  paper  is  organized  as  follows.  The  problem  formulation 
is  explained  in  Section  II.  A  backstepping-based  tracking  con¬ 
troller  for  wing  twist,  which  is  inspired  by  [30],  is  derived  in 
Section  III.  A  controller  for  wing  twist  with  actuation  at  the 
wing  tip  and  its  adaptive  version  are  presented  in  Section  IV. 
A  controller  to  track  tip  displacement  commands  is  derived  for 
the  bending  dynamics  in  Section  V,  where  an  observer-based 
approach  is  introduced.  Numerical  simulations  are  presented  in 
Section  VI.  Einally,  experimental  validation  of  the  perturbation- 
observer-based  control  of  beam  bending  is  presented  in 
Section  VII. 

II.  Motivation  and  Problem  Eormulation 

The  motivation  to  consider  the  problem  of  boundary  control 
of  beams  stems  from  the  problem  of  controlling  flexible  wings 
for  agile  aircraft  maneuvers.  Suppose  that  the  wing  has  length 
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Fig.  2.  Figure  showing  the  problems  addressed  in  this  paper. 

L,  a  mass  per  unit  span  of  m,  and  let  y  denote  the  span  wise 

coordinate.  Furthermore,  let  ^  =  y)  and  0  =  y)  denote 
the  bending  and  twist  displacements,  respectively,  as  shown  in 
Fig.  2.  Let  Ell)  and  GJ  denote  the  bending  and  torsion  stiffness, 
respectively,  where  Ii,  is  the  second  moment  of  area  of  the  cross 
section  about  the  local  bending  axis,  and  J  is  the  torsional 
constant.  Let  Ip  denote  the  polar  moment  of  inertia  of  the  wing 
cross  section. 

The  resultant  wind  velocity  at  a  given  point  on  the  wing 
u  is  a  sum  of  the  aircraft  flight  speed  (measured  in  the 
body  axis)  and  the  wing  vibration  rate  u/  =  [0  0  so  that 
u  =  Twb^b  +  u/  (^),  where  Twb  is  the  rotation  matrix  from 
the  body  frame  to  the  local  wing  frame.  Let  c  denote  the  wing 
chord  length.  Let  x^c  denote  the  distance  between  the  center 
of  mass  and  the  shear  center  of  the  wing,  and  Xa  c  the  distance 
between  the  aerodynamic  center  and  shear  center.  The  wing 
is  loaded  transversely  with  a  load  F5,  which  depends  on  the 
wind  velocity,  inertial  effects  that  arise  from  aircraft  acceleration 
and  added  mass  effect,  and  gravity.  An  aerodynamic  model  to 
calculate  F5 ,  which  is  based  on  an  ODE  approximation  for  the 
time- varying  flow  dynamics  on  the  wing,  can  be  found  in  [18]. 

Then,  the  PDE  model  of  the  structural  dynamics  of  the  right 
wing  is  given  by 


Note  that  both  root  and  tip  control  configurations  have  one  Neu¬ 
mann  and  one  Dirichlet  boundary  condition.  In  (2)  and  (3),  Ft  ip 
and  Mtip  are  the  applied  tip  shear  force  and  twisting  moment, 
respectively.  The  root  actuation  variables,  5r  and  Or,  are  the 
rigid  dihedral  (up  and  down  motion)  and  wing  incidence  an¬ 
gles.  The  term  77  denotes  the  Kelvin- Voigt  damping  coefficient. 
It  need  not  be  the  same  for  bending  and  twisting  in  general. 

The  right-hand  side  (RHS)  of  (1)  depends  on  u^,  which  is 
given  by 

m(uB +S'(u;b)ub) +m  /  (li/ +  S'(u;B)u/)dy  =  Fnet 

J  W 

Jojb  +  S{u}b)Ju}b  +  /  {Ip{y)ujf  +  S{u)B)Ip{y)ojf)dy 

J  W 

=  Mnet  (4) 

where  m  is  the  total  mass  of  the  aircraft,  J  is  the  moment  of 
inertia  tensor  for  the  aircraft,  denotes  integration  over  the 
wing,  S'(-)  denotes  the  cross  product,  and  Fnet  and  Mnet  rep¬ 
resent  the  net  external  (aerodynamic  +  gravitational)  force  and 
moment  on  the  aircraft.  Eurthermore,  ur  and  cj/  =  [0  0  0]^ 
are  vector  representations  of  the  aircraft  angular  velocity  and 
the  twist  rate  of  the  wing,  with  components  in  the  aircraft  body 
axes.  The  dynamics  of  the  motion  of  the  center  of  gravity  of 
the  aircraft,  arising  from  the  movement  of  the  wings,  have  been 
ignored.  See  [39]  for  a  detailed  and  more  general  derivation  of 
the  equations  of  motion.  The  important  point,  though,  is  that  the 
structural  dynamics  of  the  wing  are  coupled  to  the  flight  dynam¬ 
ics  of  the  complete  aircraft,  which  are  nonlinear  in  their  own 
right,  which  makes  control  design  a  challenging  assignment. 

Physically,  Ftip  and  Mtip  can  be  realized  by  using  wing  tip 
flaps,  not  unlike  the  outboard  feathers  on  the  wings  of  a  bird. 
Indeed,  trailing-edge  effectors  have  been  demonstrated  to  be 
effective  at  wing  flutter  suppression  as  well  [3].  On  the  other 
hand,  6r  and  Or  are  easier  to  control,  and  the  capability  for  root 
actuation  is  present  naturally  in  flapping  wing  aircraft. 

Problem  Formulation 

The  control  objective  is  to  ensure  that 
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where  the  subscripts  t  and  y  denote  partial  derivatives,  i.e., 
itt  =  §^,^tyyyy  =  ,  and  SO  on.  A  similar  equation  Can  be 

derived  for  the  left  wing  as  well.  The  boundary  conditions  for 
tip-based  actuation  are  given  by 

0)  =  4  {t,  0)  =  4y  {t,  L)  =0 

=  »(f,o)  =  o.  =  (2) 

Eh  GJ 

while  those  for  root-based  actuation  (also  shown  in  Pig.  1)  are 
0)  =  iyy  {t,  L)  =  0  =  ^yyy  {t,  L)  =0 
^yit,0)  =  Snit),  d{t,0)  =  enit),  eyit,L)  =  0.  (3) 


Hm  6{t,y)dy  —  H{t)^  =0  (net  lift)  (5) 

Urn  y0{t,  y)dy  —  IIi{t)^  =  0  (rolling  moment)  (6) 

lim  (  [  ^y  (t,  y)dy  -  R(t)]  =  lim  (^(t,  L)  -  R{t))  =  0 

(bending  displacement  of  the  tip)  (7) 

where  H{t),  Hi{t),  and  R{t)  are  sufficiently  smooth  time- 
varying  signals.  Although  we  state  asymptotic  convergence 
as  the  objective,  we  will  prove  exponential  convergence  to  a 
uniform  ultimate  bound.  The  net  lift  produced  by  the  wing  is 
/cf  0.5p||u||2  cGr  (t/,  0)dy,  where  p  is  the  air  density,  ||u||  is 

the  local  wind  speed,  and  Gr  is  the  coefficient  of  lift,  which  de¬ 
pends  on  the  angle  of  attack  of  the  aircraft  a,  and  the  local  wing 
twist  angle  0.  Depending  on  the  aerodynamic  theory,  Gr  may  be 
a  nonlinear  function  of  a  and  0,  and  it  is  almost  always  a  function 
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of  y  for  a  finite-span  wing  [36].  The  term  (t,  y)dy  =  ^{L) 
is  a  measure  of  the  effective  wing  dihedral,  which  is  a  key 
yaw  control  parameter  [39].  It  is  approximately  the  ratio  of  two 
components  of  .  The  first  points  are  along  the  y  direction  in 
Fig.  2,  and  the  second  points  upward.  Hence,  it  is  a  measure  of 
the  amount  of  side  {y-)  force  produced  by  the  wing  which,  in 
turn,  produces  a  yawing  moment  on  the  aircraft. 

Remark  1:  The  well-posedness  of  the  closed-loop  systems 
considered  here  can  be  shown  by  proving  that  the  input-output 
map  of  the  system  exists  and  is  bounded  [7],  [48].  For  the 
twisting  dynamics  actuated  by  root  control,  this  is  achieved  by 
designing  the  control  to  map  the  system  onto  well-posed  and 
exponentially  stable  dynamics.  For  twisting  dynamics  actuated 
at  the  wing  tip,  the  input-output  map  is  actually  a  finite-order 
ODE,  and  its  well-posedness  follows  from  the  standard  exis¬ 
tence  and  uniqueness  theorems  for  ODEs.  The  well-posedness 
of  the  closed-loop  bending  dynamics  can  be  shown  using  a 
transfer  function  approach  [7],  [12]. 

Remark  2:  It  was  shown  in  [12]  that  the  boundary  condition 
Ehiyyy  =  Ftip  in  (2)  needs  to  be  replaced  with  EIi,{^yyy  + 
Vityyy)  =  -Ftip  when  the  Kelvin- Voigt  damping  coefficient 
7/  7^  0  in  order  to  ensure  the  well-posedness  of  the  PDE.  Note, 
however,  that  this  replacement  essentially  creates  a  low-pass 
filter  for  the  control  input  Ftip.  As  long  as  the  dominant  fre¬ 
quencies  in  Ftip  are  smaller  than  the  cutoff  frequency  l/y,  the 
boundary  condition  in  (2)  is  sufficiently  accurate  for  control 
design.  Moreover,  the  Kelvin-Voigt  damping  term  is  physically 
akin  to  distributed  load  [F^  in  (1)]  on  the  wing,  which  does  not 
affect  the  boundary  conditions. 

III.  Root-Based  Tracking  Controller  for  Lift 

In  this  section,  we  consider  a  twist  PDE  with  constant  coeffi¬ 
cients.  The  PDE  backstepping  method  is  used  to  design  a  con¬ 
troller,  which  ensures  lim^^oo(/o^  ddy  —  H(t))  =0  [see  (5)]. 
The  twisting  dynamics,  with  the  angle  of  incidence  at  the  wing 
root  as  the  control  input  and  the  wing  tip  free,  are  given  by 

{t,  y)  -  bOtyy  {t,  y)  -  a0yy  (t,  y)  =  Me{t,  y) 

e{t,{))  =  u{t),  ey{t,L)  =  {).  (8) 

We  first  make  a  coordinate  transformation  and  define 

w{t,y)  =  /:  0(t,x)dx.  (9) 

Physically,  w{t,y)  measures  the  lift,  which  is  generated  by  the 
outboard  section  of  the  wing  starting  at  y  and  terminating  at 
the  tip.  Note  that  Oy  (t,  F)  =  0  at  the  free  end  y  =  L.  Hence,  it 
follows  that 


Thus,  the  dynamics  of  w  are  described  by  the  PDE 

Wtt  —  bwtyy  —  CLWyy  =  Mrc,  w{t,  L)  =  0,  Wy  (t,  0)  =  U (t) 

(11) 

where  U{t)  is  the  control  signal  [Or  in  (3)],  and  recall  that 
a  =  GJ / Ip  and  b  =  ya.  Eurthermore,  we  defined  M  so  that 
MO  =  —XacFjjj Ip,  where  F^  is  a  linear  function  of  0.  Note  that 
Wy  (t,  y)  =  0{t,  y).  The  control  objective  is  now  recast  to  ensure 

that  lim^^oo 0)  = 

The  method  presented  in  this  section  rests  on  the  following 
steps: 

1)  obtain  a  backstepping  transformation  w  ^  v,  which  maps 
the  dynamics  in  (11)  into  desired  stable  dynamics; 

2)  compute  the  boundary  condition  0);  and 

3)  derive  a  trajectory  planning-based  design  with  the  bound¬ 
ary  conditions  of  the  v  dynamics. 

Step  1  (Volterra  transform):  Define  the  Volterra  transform 

/y 

k{y,x)w{t,x)dx  (12) 

and  the  u-dynamics 

vtt—bvtyy  —  avyy  =  Mv  —  bpvf  —  apv^  u(f,  F)  =  0,  p  >  0. 

(13) 

Note  that  we  have  prescribed  only  one  boundary  condition  for  v, 
which  follows  from  the  fact  that  w{t,L)  =0.  The  second  bound¬ 
ary  condition  u^(t,0)  has  to  be  calculated  to  ensure  that  the 
lim^^oo  ^(F  0)  ~  ~H  (t)  (this  is  what  we  can  guarantee  in  prac¬ 
tice,  in  place  of  asymptotic  convergence).  The  target  dynamics  in 
( 1 3)  are  very  similar  to  those  designed  for  the  regulation  problem 
in  [30,  ch.  7,  pp.  79-89].  However,  the  boundary  conditions  are 
chosen  to  match  the  control  objective  in  (5).  Eor  the  regulation 
problem  (H  =  0),  we  can  prescribe  Vy  (f,  0)  =  v{t,L)  =  0  [30]. 
In  contrast,  in  order  to  track  a  nonzero,  time-varying  i7(t),  we 
choose  Vy  (F  0)  by  using  trajectory  planning. 

We  first  establish  the  stability  of  (13)  with  homogeneous 
boundary  conditions,  i.e.,  Vy{t,  0)  =  v{t,  L)  =  0. 

Proposition  1:  The  target  dynamics  (13)  with  u^(FO)  = 
v{t,L)  =0  are  stable  if  and  only  if 

f  M  TT^  TT^  \ 

(14) 

Proof:  Using  the  method  of  separation  of  variables,  it  is  easy 
to  check  that  the  eigenvalues  of  the  resulting  system  are  the 
solutions  of 


ry 

tit,y)  =  J 

/: 


u)tt[t,y)  =  j  9tt(t,x)dx 


{bOtxx  {t,  x)  +  aOxx  {t,  x)  +  M0{t,  x))dx 


=  b0ty{t,y)  +  a0y{t,y)  +  Mw{t,y) 

=  bwtyy{t,y)  +awyy{t,y)  +  Mw{t,y). 


A,^  +  b{i'^  +  p)X  +  +p)  —  M)  =  0 

,  2n  +  1  TT 

where  u  = - — ,  n  =  0,1,2, _ 

Z  1j 


(15) 


(10) 


The  target  dynamics  are  stable  if  and  only  if  the  control  design 
parameter  p  satisfies  (14).  ■ 

Next,  we  solve  iork{y^x)  in  the  Volterra  transform  (12). 


PARANJAPE  et  al.\  PDE  BOUNDARY  CONTROL  EOR  ELEXIBLE  ARTICULATED  WINGS  ON  A  ROBOTIC  AIRCRAET 


629 


Proposition  2:  The  gain  kernel  k{y ^x)m{\2)  is  the  solution 
of  the  Klein-Gordon  PDE 


kxx{y,x)  -  kyy{y,x)  =  -pk{y,x),  k{y,L)  =  0 
k{y,y)  =  ^{L-y)  (16) 

and  given  in  terms  of  the  modified  Bessel  function,  Iii'Y 


Proof:  We  can  derive  expressions  for  Vu  and  Vyy  from  (12) 

vtt  =  wtt  -  J  k{y,x)wtt{x)dx 


-  bWfyy  T  CLWyy  T  ]\d  W 

ry 


ry 

J  k{y,  x)  {bwtxx  {x)  +  awxx  {x)  +  Mw{x))  dx 


-  bWfyy  T  CLWyy  T  ]\/[  W  - 

y 


/; 


k{y^  x)Mw{x)dx 


I 


kxx  (2/,  x)  {bwt  (x)  +  aw{x))  dx 

-  k{y,y){bwty{y)  +  awy  (y))  +  k^  {y,y){bwt{y)  +  aw{y)) 
+  k{y,L){bwty{L)  +  aWy{L)).  (18) 

The  expression  for  Vyy  can  be  derived  as  follows: 

ry 

Vy  [y)  =  Wy  [y)  -  k{y,  y)w{y)  -  ky  {y,  x)w{t,  x)dx 
Vyy  {y)  =  Wyy  {y)  -  2ky  (j/,  y)w (y) 

-  kx(y,y)w{y)-k(y,y)wy{y)- j  kyy{y,x)w{t,x)dx. 

(19) 


The  expression  for  Vtyy  is  similar  to  that  for  Vyy.  We  substitute 
(18)  and  (19)  into  (13)  and  isolate  the  coefficients  of  w  and 
Wy  in  the  integrand  as  well  as  outside  the  integral.  This  yields 
(16).  ■ 
Step  2  (Boundary  condition,  u(t,  0),/<9r  tracking):  We  have 
already  ascertained  that  i;(t,  L)  =  0.  From  (12),  it  follows  that 

u(t,  0)  =  u;(f,  0)  +  f  k{0^x)w{t,x)dx  (20) 

Jo 

which  does  not  yield  an  exact  relationship  between  H (t)  de¬ 
fined  in  (5)  and  u(t,  0).  Instead,  we  need  to  determine  an  ap¬ 
propriate  reference  value  for  denoted  by  Hr(t),  which 

ensures  that  w{t,0)  approximately  tracks  —H(t).  For  exam¬ 
ple,  we  could  approximate  w{t,  x)  fi  {x)w{t,  0),  fi  (0)  =  1, 
where  fi  (x)  denotes  the  shape  of  the  first  twisting  mode.  Thus, 
weseti^^(t)  =  -h  Jq  /c(0,  x)/i  (x)(ix).  It  is  important 

to  appreciate  that  ignoring  the  higher  modes  has  no  repurcus- 
sions  for  stability  since  it  is  guaranteed  separately.  At  the  same 
time,  it  may  be  necessary  to  include  more  than  one  mode  if  there 


modified  Bessel  function  In(y)  satisfies  y‘^ +  yln(y)  —  (y‘^  + 
n‘^)In(y)  =  0. 


are  tight  bounds  on  the  tracking  error.  A  similar  approach,  which 
is  based  on  backstepping,  but  for  constant  reference  signals,  has 
been  derived  in  [49] . 

Step  3  (Trajectory  planning  for  wing  twist):  We  design  a  tra¬ 
jectory  planning-based  algorithm  (see  [30],  Ch.  12).  We  define 
a  new  state  ,  where  the  superscript  “r”  denotes  the  reference 
value.  The  dynamics  of  are  given  by  a  PDE  that  is  identical 
to  (13): 

-  Kyy  -  o-vly  =  Mv^  -  bpvl  -  apv^ ,  {t,  L)  =  0 

(21) 

where  Vy  (t,  0)  is  the  control  input  to  be  designed  using  trajectory 
planning  to  ensure  that  (t,  0)  tracks  Hp  (t). 

Since  (f,  L)  =  0,  we  expand  (f,  y)  as  a  polynomial 

v''  {t,  y)  =  y2  0)  -f  ■  (22) 

i=i 

Substituting  for  into  (21),  we  get 

VjO)  +  bPVjiO  +  {ap  -  M)rij{t)  =  677^+2 (i)  +  arij+2{t). 

(23) 

The  requirement  that  v’’  {t,  0)  =  Hr  (t)  gives 

^  Li 

=  Hr{t).  (24) 

i=i 

The  value  of  N  can  be  chosen  on  a  case-by-case  basis.  As  an 
illustration,  if  we  truncate  the  series  at  =  3,  we  get 

Vi  {t)  +  bpfii  (t)  +  {ap  -  M)rii  (t)  =  bris{t)  +  arisit).  (25) 

A  similar  equation  can  be  obtained  for  772  (f),  with  the  RHS  equal 
to  zero.  Therefore,  p2  converges  to  zero  exponentially  fast  and 
can  be  ignored.  The  requirement  that  (t,  0)  =  Hp  (t)  gives 

111  {t)L  +  773  (t)  ^  =  Hr  {t) .  (26) 

6 

Substituting  7^3 (t)  =  ^  {Hp{t)  —  into  (25)  gives 

771  (t)  +  b(p+  771  (t)  +  ^ap-M+  771  (t) 

=  ^  (^bHriO+aHriOf  (27) 

Then,  Vy  (t,  0)  is  calculated  by  differentiating  both  sides  of  (22) 

{t,0)  =  -111  (t)  -  773  (t)  k  =  2771  (f)  -  (28) 

where  vl{t,y)  =  dv^{t,y)/dy. 

Proposition  3:  Based  on  Steps  1-3,  we  now  define  U{f)  = 
6>(t,  0)  =  Wy{t^  0)  as  illustrated  in  the  block  diagram  in  Fig.  3. 
First,  we  set  the  boundary  condition  7;^(t,0)  =  7;J(t,0).  Let 
v{t^y)  =  v{t,y)  —  (t^y).  It  is  straightforward  to  check  that 

the  dynamics  of  v  are  given  by  the  exponentially  stable  PDE 

Vtt  -  bvtyy  -  avyy  =  Mv  -  bpvt  -  apv 

v{t,L)  =  Vy{t,Q)  =  0.  (29) 

Thus,  in  particular,  it  follows  that  v{t,  0)  ^  (f ,  0)  exponen¬ 

tially  fast  and  0)  tracks  Hp  (t). 


630 


IEEE  TRANSACTIONS  ON  ROBOTICS,  VOL.  29,  NO.  3,  JUNE  2013 


Vy^m 


V^(t,0) 


w  dynamics 


w(t,0) 


w(t,y) 


v''  dynamics 


■H(t) 


Fig.  3.  Block  diagram  showing  the  tracking  controller  for  twist  (6)  dynamics. 
This  structure  is  identical  to  the  classic  strict  feedback  structure  for  systems 
described  by  ODEs  [29]. 


An  expression  for  Vy  (t,  y)  can  be  found  by  by  differentiating 
both  sides  of  the  Volterra  tranform  in  (12) 

ry 

Vy  (t,  y)  =  Wy  (t,  y)  -  k{y,  y)w{t,  y)  -  ky  {y,  x)w{t,  x)dx. 

(30) 

Finally,  the  control  input  f/(f)  is  calculated 

U  (t)  =  0(t,  0)  =  Wy  (f,  0) 

=  u^(f,  0)  + /c(0, 0)u;(t,  0)  —  f  ky{0,  x)w{t,  x)dx  (31) 

Jo 

^ - V - " 

feedback 

where  Vy  (t,  0)  ^  Vy  (t,  0)  exponentially  fast.  The  term  Vy  (t,  0) 
is  a  feed- forward  tracking  term  in  the  control  signal,  while  the 
remaining  two  terms  act  as  feedback  stabilizers. 

We  now  assert  that  the  closed-loop  system  is  stable  in  the 
sense  of  £2.  Given  a  quantity  y)  G  M,  define 

IMOWc,  =  ]J  VL{t,yfdy. 

Theorem  1:  The  closed-loop  rc-dynamics  (11)  and  (31)  are 
stable  in  the  sense  of  £2,  i.e.,  ||fc(t)  ||/:2  is  bounded. 

Proof:  Using  a  Lyapunov  stability  argument,  it  can  be  shown 
that  ||f^(t)||/:2  is  bounded  at  all  times.  Since  (12)  is  a  dif- 
feomorphism,  we  can  find  another  function  l{y^x)  (which  is 
expressed  in  terms  of  the  Bessel  function  Ji  [30])  such  that 
w{t,  y)  =  v{t,  y)  —  Kv^  x)v{t,  x)dx,  and  l{y,  x)  is  bounded 
for  alio  <y<x<  L.  Since  the  ||f^(f)||/:2  is  bounded,  it  follows 
that  II  u;(t)  1 1  £2  is  bounded.  ■ 

Remark  3:  The  approach  that  is  detailed  in  this  section  also 
works  when  the  RHS  is  nonlinear,  i.e.,  M  =  M{y).  Lineariza¬ 
tion  can  be  used  in  the  general  case,  where  the  RHS  is  of  the 
form  M(^, 6>i, 6^2).  It  has  been  shown  in  [38]  that  a  regulator 
based  on  (3 1)  can  stabilize  the  twisting  dynamics  in  the  presence 
of  these  nonlinearities.  The  analysis  is  not  repeated  here. 
Remark  4:  A  few  observations  are  worth  noting  here. 

1 )  For  stability,  it  is  essential  that  b  >  0,  i.e. ,  the  Kelvin-Voigt 
damping  coefficient  is  always  positive.  The  aerodynamics 
could  introduce  negative  damping,  but  it  can  be  compen¬ 
sated  by  the  term  bpvf ,  and  wing  flutter  can  be  prevented. 


The  need  to  compensate  for  negative  damping  imposes  an 
additional  constraint  on  p. 

2)  Damping  and  stiffness  cannot  be  added  independently. 
They  are  added  in  the  ratio  b/a. 

3)  The  controller  in  (31)  requires  that  the  twist  angle  at  all 
points  on  the  wing  be  known.  This  difficulty  can  be  cir¬ 
cumvented  by  designing  a  PDF-based  observer  [30]  or, 
practically,  using  a  series  of  distributed  sensors  and  fitting 
their  output  with  an  a  priori  designed  spline. 

4)  Finally,  the  control  law  does  not  require  that  a,  b,  or  M 
be  known  for  the  purpose  of  regulation.  We  only  need  to 
know  bounds  on  a  and  M  to  choose  the  gain  p. 

IV.  Boundary  Control  of  Twisting  Motion: 

Wing  Tip  Control 

Although  backstepping  can  be  employed  when  the  twisting 
moment  at  the  wing  tip  is  used  as  the  control  input,  it  turns  out 
that  there  is  a  simpler  alternative,  as  described  in  this  section.  It 
relies  on  the  fact  that  the  system  has  a  finite  relative  degree  for  the 
input-output  combination  of  the  tip  moment  and  y)dy, 

respectively. 

A.  Tip  Boundary  Control  When  All  Parameters  Are  Known 

As  in  the  previous  section,  one  can  design  a  backstepping 
controller  for  the  case  where  a  control  moment  is  applied  to 
the  free  end  {y  =  L)  of  the  wing,  while  the  other  end  {y  =  0) 
is  clamped.  In  fact,  the  procedure  in  both  cases  is  identical,  al¬ 
though  the  final  expressions  for  the  control  law  differ  slightly. 
Alternately,  in  the  case  of  MAVs,  one  may  do  without  a  sta¬ 
bilizing  controller.  The  “tracking  half”  of  the  controller  may 
be  designed  by  using  the  output  measurements.  This  method  is 
useful  for  adaptive  designs  as  well.  We  consider  the  wing  model 

On  -  bOtyy  -  aOyy  =  M9,  9{t,  0)  =  0,  9y  {t,  L)  =  u{t) 

(32) 

where  the  control  input  is  a  moment  that  is  applied  at  the  wing 
tip,  {u{f)  =  0y{L,f)  =  Mtip  in  (2),  and  b  =  rja  =  r]GJ /Ip. 
Furthermore,  we  defined  M  so  that  MO  =  —XacFi)/ Ip,  where 
F5  is  a  linear  function  of  0. 

Theorem  2:  A  dynamic  controller  of  the  form 

bu{f)  +  au{f)  =  H{f)  -  MH{t)  -  (M  -h  k)e{t)  -  kce{t) 

+  a0y{f0)Fb0ty{f0)  (33) 
ensures  the  following  control  objective: 

^lim  9{t,y)dy  -  H{t/j  =0.  (34) 

Proof:  Let  e{t)  =  0{t,y)dy  —  H{t)  denote  the  error, 
which  needs  to  be  regulated.  Then, 

e  =  [  0u{t,y)dy  -  H{t) 

JO 

L 

{ciOyy  +  bOtyy  +  MO)  dy  —  H{t) 

=  aOy  (t,  L)  -  aOy  (t,  0)  -h  bOty  (t,  L)  -  bOty  (t,  0) 
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+  Me(t)  -  H{t)  +  MH{t) 

=  hu{t)  +  au{t)  —  aOy  (t,  0)  —  bOfy  (t,  0) 

+  Me{t)  -  H{t)  +  (35) 

The  controller  (33)  renders  the  system  into  a  globally  expo¬ 
nentially  stable  linear  system  with  kc  >  0,  k  >  0 

e{t)  +  kccit)  -h  ke{t)  =  0.  (36) 

Hence  proved.  ■ 

Remark  5:  The  control  law  in  (33)  suggests  that  0  need  not  be 
monitored  or  measured  at  all  locations  on  the  wing.  Instead,  only 
Oy  (t,  0)  needs  to  be  measured  or  estimated.  The  reference  signal 
H (t)  is  known.  It  may  be  difficult  to  inject  damping  because  e{t) 
is  the  rate  of  change  of  the  lift  and,  in  practice,  would  require 
differentiating  noisy  acceleration  signals.  Another  interesting 
observation  is  that  although  the  PDE  system  had  an  infinite 
relative  degree  when  the  root  twist  was  chosen  as  the  control 
input,  the  relative  degree  is  2  when  twisting  moment  at  the  wing 
tip  is  considered  as  the  input.  Finally,  the  control  law  design  that 
is  described  is  this  section  lends  itself  readily  to  adaptation  if  a 
and/or  M  are  unknown. 


by  M/a  <  jL^ ,  it  follows  that  the  stability  margin 

improves  by  a  factor  of  16  using  only  the  tracking  controller, 
although  it  does  not  stabilize  the  wing  for  all  values  of  M  and 
a  as  backstepping  does.  In  principle,  the  tracking  controller 
converts  the  wing  from  a  cantilever  beam  to  a  clamped-clamped 
beam.  In  practice,  this  translates  to  the  ability  to  increase  the 
wing  flexibility  by  an  order  of  magnitude,  or  increase  the  wing 
divergence  speed  fourfold. 

C.  Adaptive  Control  for  a  Wing  Tip  Actuator 

Consider  the  dynamics  in  (32)  with  the  objective  in  (34),  and 
suppose  that  M  is  unknown.  The  control  law  in  (33)  is  modified 
so  that 

bu{t)  +  au{t)  =  H(t)  —  M(t){H(t)  +  e(t))  —  ke{f) 

-  kce{t)  +  a0y{t,0)  +  b0ty{t,0)  (41) 

where  M{t)  is  the  estimated  value  of  M,  and  e(t)  = 
0{t^y)dy  —  H (t)  denotes  the  tracking  error.  The  closed-loop 
tracking  error  error  dynamics  are  described  by  the  ODE 


B.  Tracking  and  Stability 

The  problems  of  tracking  and  stabilization  are  distinct  be¬ 
cause  the  PDE  system  (32)  is  infinite  dimensional.  Nevertheless, 
a  tracking  controller  improves  stability  as  described  presently. 
One  way  to  measure  the  influence  of  a  tracking  controller  on  the 
stability  of  the  closed-loop  system  is  to  set  iif  =  0.  Make  a  co¬ 
ordinate  transformation  w{t,  y)  =  0{t,  z)dz,  similar  to  (9), 

so  that  achieving  =  0  is  equivalent  to  achieving  u;(t,  0)  =0. 
It  follows  that  the  transformed  dynamics,  with  H  =  ^  obtained 
by  the  tracking  control  algorithm,  are  given  by 

Wtt  (c  y)  -  bwtyy  (t,  y)  -  aWyy  (t,  y) 

=  Mw{t,  y)  -  bwtyy  (t,  0)  -  awtyy  (t,  0) 

w{t^  L)  =  w{t^0)  =  0,  w'{t,0)=0.  (37) 

The  third  boundary  condition  is  not  entirely  independent.  Let 
w{t,  y)  =  irj(t)(j){y).  Then,  we  get 

fi{t)  -  Mrjjt)  ^  (t)yy{y)  -(pyyjO)  ^  _^2 
bf]{t)  +  aT]{t)  (f>{y) 

where  k  is  a  constant.  It  follows  that  the  condition  for  stability 
(exponential  stability  to  be  precise,  since  the  present  system  is 
linear)  is  M/a  <  k‘^.  The  differential  equation  for  ^{y)  in  (38) 
can  be  solved  to  get 

(p{y)  =As'm{ky)  +  Bcos{Xy)  +  4>yy{0)/k‘^.  (39) 

The  boundary  conditions  0(0)  =  0^(0)  =  0  lead  to  0(^)  = 
B{cos{ky)  —  1).  Finally, 

0(1/)  =  0  k  =  2n7r/L^  0(1/)  =  0,  n  =  l,2, - 

(40) 

Had  we  not  imposed  the  condition  77(0)  =  0,  we  would  have 
obtained  M/a  <  7r^/(4I/^)  as  the  condition  for  (exponential) 
stability.  Since  the  condition  for  closed-loop  stability  is  given 


e{f)  -h  kce{f)  -h  ke(f)  =  —M{t){H{t)  +  e(t))  (42) 

where  M(t)  =  ]M(t)  —  M{t).  We  design  a  projection-based 
adaptive  law  for  M(t)  to  ensure  that  the  error  e{t)  remains 
bounded  with  as  small  a  bound  as  possible.  The  projection  op¬ 
erator,  which  is  denoted  by  Proj(',  •),  is  used  commonly  in  the 
adaptive  control  literature  [6],  [35]  for  parameter  estimation. 

Theorem  3:  Let  x{f)  =  [e(t),  e(t)]^ .  Suppose  that  there  exist 
positive  constants  Bi  and  B2  such  that  \M\  <  Bi  and  \M\  < 
B2 .  Then,  the  following  adaptive  law 


M{t)  =  7Proj  (  M{t),x^P 


{e{t)  +  H) 


(43) 


with  \M{t)\  <  Bi  ensures  that  there  exists  a  positive  constant 
K {Bi ,  B2 )  and  a  time  T  >  0  such  that 

||x(i) lloo  <  \\x{t)\\2  <  t  >  T.  (44) 

\/7 


Proof:  see  [38]. 

Remark  6:  The  bound  on  ||x(t)||oo  can  be  made  arbitrarily 
small  by  choosing  a  large  7.  The  steady-state  beam  shape  of 
the  wing  depends  on  the  steady- state  value  of  the  error  M{f). 
Finally,  it  is  worth  noting  that  although  a  and  b  were  assumed 
to  be  known,  the  aforementioned  analysis  can  be  repeated  to 
accommodate  an  unknown  a  and  b  as  well. 


D.  Control  of  Rolling  Moment 

An  abstract  measure  of  the  rolling  moment  is  y0{t^  y)dy, 
which  is  defined  in  (6).  Let  ept^y)  =  y0{t^y)dy  —  Hpt), 
where  Hi  (t)  denotes  the  reference  value  of  the  rolling  moment 
to  be  tracked,  and  ei  (t)  is  the  tracking  error.  Differentiate  ei  (t) 
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twice  with  respect  to  time 

ei(0=  [  yffu{t,y)dy  -  Hi{t) 

Jo 

L 

y{bOtyy  {t,  y)  +  aOyy  {t,  y))dy 

+  Mf  ye{t,y)dy-Hi(t) 

Jo 

=  LibOty  {t,  L)  +  aOy  {t,  L))  -  b{0t  {t,  L)  -  8t  (t,  0)) 

-  a{e{t,  L)  -  e(t,  0))  +  Mei{t)  + 

(45) 

An  interesting  observation  is  that  ei  (t)  has  a  relative  degree  of 
2  with  respect  to  0y{t,L)  (tip  control)  as  well  as  6>(0)  (root 
control).  Therefore,  a  considerably  simpler  controller  than  the 
backstepping  controller,  on  the  lines  of  the  adaptive  controller 
in  the  previous  section,  can  be  implemented  for  a  root-based 
actuator  as  well.  Indeed,  in  aicraft,  the  lift  is  controlled  using 
the  horizontal  tail,  while  wing-based  flaps  (ailerons  and  spoilers) 
are  used  primarily  for  roll  control. 

V.  Root  Boundary  Control  of  Wing  Bending 

In  this  section,  we  first  determine  the  conditions  under  which 
there  is  a  significant  time- scale  separation  between  the  bending 
and  twisting  dynamics.  When  the  time-scale  separation  is  suf¬ 
ficient,  the  control  laws  for  bending  and  twist  can  be  designed 
separately,  which  simplifies  the  design  process  considerably  in 
practice.  Time-scale  separation  is  used  quite  routinely  to  design 
flight  control  laws  [24],  [45],  [53],  [59].  We  introduce  one  of  the 
main  results  of  this  paper — the  new  perturbation-observer-based 
method  to  design  a  bending  controller.  In  particular,  we  design 
a  root-based  bending  controller  [6r  in  (3)]  using  a  combina¬ 
tion  of  the  perturbation-observer-based  approach  and  trajectory 
planning.  The  design  of  a  tip  controller  [Ft ip  in  (2)]  is  straight¬ 
forward  using  the  approach  that  is  established  in  this  section, 
and  we  demonstrate  it  experimentally. 

The  bending  PDE  in  (1)  can  be  written  as 

V^^b^tyyyy  ^^b^yyyy  —  •>  ^  ^tt) 


approximate  relation,  when  the  speed  is  less  than  the  torsional 
divergence  speed 


uj0  1.5  L 
1  -|-  u  c 


(47) 


where  uoe  and  denote  the  angular  frequencies  of  torsion  and 
bending,  respectively,  u  denotes  the  Poisson  ratio  (typically  in 
the  range  0. 1-0.4),  and  the  ratio  L/cis  called  the  aspect  ratio  of 
the  wing  (of  the  half  wing,  to  be  precise).  Therefore,  the  twisting 
dynamics  are  faster  than  the  bending  dynamics,  and  moreover, 
the  time- scale  separation  increases  with  the  aspect  ratio,  and 
reduces  with  increasing  flight  speed  [39].  In  most  practical  sit¬ 
uations,  a  time- scale  separation  of  3-5  is  considered  sufficient 
to  design  and  implement  decoupled  control  laws,  i.e.,  decou¬ 
pled  controllers  for  bending  and  twist  can  be  safely  designed 
for  wings  with  aspect  ratio  Ljc  >  3. 

Remark  7:  Wings  with  smaller  aspect  ratios  would  have  to 
be  modeled  as  plates  rather  than  beams.  Although  several  MAV 
wings  do  have  aspect  ratios  smaller  than  3  (such  as  the  one  built 
by  the  authors  [40]),  their  structure  consists  of  a  skin  wrapped 
around  or  on  a  more  traditional  beam-like  structure.  The  control 
design  explained  in  this  paper  can  still  be  employed  for  such 
wings.  The  two-time-scale  approach  is  rigorously  applicable  to 
coupled  infinite-dimensional  systems,  as  illustrated  for  thermal- 
fluid  dynamics  in  [57],  and  for  flexible  structure  dynamics  in 
[51].  Moreover,  although  the  coupled  dynamics  overlap  over  an 
infinitely  wide  frequency  range,  the  slower  subsystem  has  only  a 
tiny  fraction  of  its  energy  in  the  (semi-infinite)  frequency  range, 
where  it  overlaps  with  the  faster  dynamics  (whose  entire  energy 
is  in  that  range). 


B.  Open-Loop  Stability  of  Bending 

Assume  that  the  boundary  conditions  are  homogeneous,  i.e., 
^(t,0)  =  iyit.tS)  =  iyy{t,L)  =  ^yyy{t,L)  =  0  for  all  t.  Let 
Ob  =  EI}j/m  and  65  =  pa^.  Let  F{t,y)  =  Fn{t,y) /m,  where 
Fn  (t,  y)  succinctly  denotes  the  RHS  of  (46).  Assume  that  the 
flight  speed  Hu^  ||  is  constant. 

Proposition  4:  Subject  to  the  aforementioned  assumptions, 
the  bending  dynamics  (46)  with  =  0  are  stable  if 


^(t,0)  =  4^(t,F)  =  4^^(t,F)  =  0,  4(t,0)='u(t)  (46) 

where  the  control  input  u{f)  =  5r  (t)  in  (3)  is  designed  to  ensure 
that  i{fL)  =  R{t)  mil).  The  acceleration  term  corresponding 
mXecO  in  (1)  has  been  moved  to  the  RHS  and  merged  into  F5 
so  that  F^  =  F5  +  mXecO.  The  RHS  (46)  is  independent  of 
and  0  is  an  average  value  that  is  obtained  from  the  faster  twist 
dynamics.  Therefore,  unlike  the  twisting  dynamics,  the  onset 
of  instability  in  the  bending  dynamics  will  correspond  to  the 
damping  becoming  negative. 


Proof:  Consider  the  Lyapunov  function  (in  this  case,  the  total 
£2  energy  of  the  beam) 

v{t) = vx  y)  =  \[  dy-  (48) 


A.  Time-Scale  Separation  in  Microaerial  Vehicles  With 
Flexible  Wings 

In  [39],  it  was  shown  that  the  time-scale  separation  between 
the  twisting  and  bending  dynamics  is  given  by  the  following 


It  is  quite  straightforward  to  derive  the  following  expressions 
forF(t): 

V{i)  =  -b,  f  ilydy  +  f  y)dy.  (49) 
7o  7o 
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A  repeated  application  of  Poincare’s  inequality^  yields 

V{t)  iUy  +  £  mt, y)dy. 

Using  the  Cauchy-Schwarz  inequality,  we  get 


then  the  bending  dynamics  are  stable.  ■ 

The  expression  for  F{t,y)  in  the  simplest  case  is  of  the  form 
F{t^y)  =  where  >0  are  essen¬ 

tially  constant.  If  Cl  >  0  (a  reasonable  assumption),  it  follows 
that  niCiit  injects  damping  into  the  system  and  counters  the  po¬ 
tentially  destabilizing  contribution  from  ^2^1^  for  sufficiently 
small  values  of  However,  if  Cl{0)  is  not  uniformly  posi¬ 
tive,  and  if  it  is  out-of-phase  with  F(f,  y)  could  potentially 
cause  the  damping  to  become  negative,  leading  to  oscillatory 
instability  in  the  bending  dynamics  and  trigger  wing  flutter. 

Furthermore,  since  ki  oc  ||ul||,  the  flight  speed,  it  follows 
that  (50)  can  be  used  to  derive  a  stable  envelope  for  the  flight 
speed  and  Young’s  modulus.  However,  unlike  the  analysis  for 
wing  twist  in  Section  V-A  and  [39],  the  analysis  is  considerably 
more  complicated  because  of  the  potentially  strong  nonlinear 
contribution  from  at  low  flight  speeds. 

C.  Perturbation  Observer  for  Root  Control  of  Bending 

In  this  section,  we  design  a  high-gain  observer-based  con¬ 
troller  to  facilitate  a  trajectory-planning-based  tracking  con¬ 
troller  for  bending.  The  perturbation  observer  does  not  predict 
the  system  states.  It  uses  projection-based  adaptation  to  esti¬ 
mate  F{t,y)  which  would  be  unknown  in  practical  situations, 
although  one  may  estimate  its  spatial  profile  from  the  wing 
geometry  [36].  The  perturbation  observer  is  split  into  a  “partic¬ 
ular”  component  and  a  “homogeneous”  component  (the  notions 
will  be  made  more  precise  in  this  section).  In  particular,  the 
homogeneous  component  is  stable  and  not  driven  directly  by 
external  feedback.  Thus,  it  is  simpler  to  design  a  “tracking” 
control  law  for  it.  The  same  control  signal  is  sent  to  the  actual 
system,  whose  states  then  converge  exponentially  to  a  bounded 
envelope  around  the  observer  states.  Fig.  4  shows  a  block  dia¬ 
gram  of  the  perturbation  observer-based  controller.  It  is  the  PDE 
analog  of  feed-forward  and  feedback  control  for  systems  that 
are  described  by  ODEs.  Notice  that  the  reference  input,  R{t) 
from  (7),  enters  only  the  homogeneous  component  of  the  per¬ 
turbation  observer,  while  the  feedback  from  the  actual  system 
only  enters  the  particular  component. 

^Poincare’s  inequality  w‘^dx  <  {0)  +  w‘ldx. 


Fig.  4.  Block  diagram  for  the  perturbation  observer  coupled  to  the  system 
dynamics.  The  control  signal  u{t)  is  generated  from  trajectory  planning,  while 
R(t)  is  the  desired  reference  signal  from  (7). 

Let  F(t,  y)  =  W{t)^  fbiv)  +  cr(t),  where  IU(f)  andcr(t)  are 
unknown  and  bounded  with  known  bounds.  We  assume  that  W 
and  &  are  also  bounded  with  known  bounds.  The  set  of  functions 
(j)}j  {y)  can  be  chosen  to  get  a  satisfactory  bound  on  a,  and  using  a 
knowledge  of  the  wing  geometry  [36]  (e.g.,  the  force  distribution 
on  an  elliptical  wing  is  rectangular,  and  vice  versa). 

The  perturbation  observer  for  the  bending  dynamics  in  (46)  is 
designed  as  a  combination  of  two  components,  which  are  called 
the  “particular”  component  and  the  “homogeneous”  component, 
for  reasons  that  will  become  apparent  presently.  Let  and 
denote  the  states  of  the  particular  and  homogeneous  compo¬ 
nents,  respectively.  The  dynamics  of  ^p  and  are  described  by 
the  following  PDEs: 

Cpdt  F  ^h^p^tyyyy  ifi  V)  F  CLh^p^yyyy  ifi  V) 

=  -bbpip,t{t,y)  -  abpip{t,y)  +W{t)^(pb{y)  +  d{t) 

ip^yy{t,L)  =  ip,yyy{t,L)  =  =  ip,y{t,0)  =  0  (51) 

for  the  particular  component,  where  Cp  =  tp  ~  and 

^h^tt  3“  ^b^h^tyyyy  3“  ^b^h^yyyy  ^bP^h^t  ^bP^h 
^h,yy{t^L)  =  ^h^yyyit^  F)  =  <^/i(f,0)  =  0,  ^}i,y(fiF)  =  u{f) 

(52) 

for  the  homogeneous  component.  Note  that  the  homogeneous 
half  is  a  stable  system,  and  its  dynamical  equation  does  not 
depend  on  ^  or  ^p  (which  is  why  it  is  called  the  “homogeneous” 
component).  Moreover,  the  reference  signal  R{f)  is  sent  only  to 
the  homogeneous  component. 

We  now  put  together  the  two  components.  Define  ^  =  tp  F 

•  Then,  the  dynamics  of  ^  are  given  by 

y)  ^bityyyy  y)  F  db^yyyy  y) 

=  (pbiy)  +  d{t)  -  bbplt{t,y)  -abpi{t,y) 

iyy{t,L)  =  iyyy{t,L)  =  i{t,0)  =0,  iy{t,0)=U{t)  (53) 

where  ^  “  6  and  p  >  0  is  chosen  to  ensure  desirable 

convergence  properties.  Recall  that  a&  =  ER/m  and  65  = 
r]ER/m.  Choose  projection-based  adaptive  laws  for  W (t)  and 
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a{t) 

W{t)  =  7„Proj  -  j\it  +  S0My)dy 

a{t)  =  7aProj  (It  +  Si)dy^  (54) 

where  7^  >  0  is  the  adaptation  rate,  and  the  bounds  on  W  and 
(j  are  chosen  to  be  equal  to  the  known  bounds  on  W  and  cr, 
respectively. 

We  are  now  ready  to  state  our  main  result,  which  asserts  that 

which  is  the  sum  of  the  states  of  the  two  observer  components, 
tracks  i.e.,  the  system  state. 

Theorem  4:  The  error  between  the  system  (46)  and  the  per¬ 
turbation  observer  (53),  in  the  sense  of  ||^(t,  y)\\c2  ^  is  globally 
uniformly  bounded,  and  the  bound  can  be  made  arbitrarily  small 
by  increasing  7^ . 

Proof:  To  streamline  the  proof,  we  introduce  some  useful 
notation.  We  define  numbers  5i ,  52 ,  ,  ^2 ,  which  depend  on  the 

gain  p,  as  follows: 


51  =  min  {{bbp  -  5),  SabP,  Sab} 

52  =  imax{l  +  6,  {ub  +  Sbb),  {{ab  +  5bb)p  +  5)} 

TP^TP  +  cr2<^,  |TP^tT|  +  |c7(t|  <  y  (55) 

where  ^  >  0  is  an  arbitrarily  small  constant,  which  is  introduced 
in  the  proof.  Equation  (55)  reflects  the  fact  that  W (t),  cr(t),  their 
time  derivatives,  and  predicted  values  are  bounded  with  known 
bounds. 

Consider  the  Lyapunov  function 

v{t)  =  \j^  (|2  +  (a,  +  SbbVeyy  +  {ab  +  5bb)p'e)  dy 

+  T  (w^W  +  <72)  +  T  iibdy  (56) 

7a  ^  ^  Jo 


where  W  =  W  —  W,  and  a  =  a  —  a.  The  constant  5  is  chosen 
to  be  small  enough  so  that  the  Lyapunov  function  is  positive 
definite.  Differentiating  both  sides  with  respect  to  time  and  per¬ 
forming  integration  by  parts  a  few  times,  we  get 


V{t)  <  -{bbP  -  S)  lldy  -  SabP  t  'edy 
JO  Jo 


pL  pL 

-Sab  /  llydy-hb  /  Ityy 
Jo  Jo 


dy 


-(w{tfW{t)  +  a{t)&{t)).  (57) 


la 

Define  a  constant  Vb  =  —  +  Lrom  (57),  we  get 

^  la  Sl7a  V  to 


V{t)  <  --P(i)  +  7  +  fill 

52  la  S2la 


■^(P(i)-Po).  (58) 

52 


It  follows  that  V (t)  converges  to  Vb  exponentially  fast.  Thus, 
the  error  dynamics  between  the  perturbation  observer  and  the 


actual  system  are  globally  uniformly  bounded,  and  the  bound 
reduces  uniformly  with  increasing  7^ .  ■ 

The  two  halves  of  the  observer,  put  together,  thus  function 
as  a  standard  perturbation  observer  whose  job  is  to  estimate  the 
unknown  nonlinearities  and  disturbances  using  full  state  feed¬ 
back.  In  standard  adaptive  control  methodologies,  a  control  law 
is  first  designed  assuming  a  complete  knowledge  of  the  system. 
Thereafter,  the  unknown  system  parameters,  which  appear  in  the 
control  law,  are  replaced  by  estimated  values  that  are  supplied 
by  the  adaptation  law. 

In  contrast,  we  decompose  the  original  system  into  a  “par¬ 
ticular”  and  a  “homogeneous”  half,  whose  states  are  predicted 
by  the  two  components  that  are  described  briefly;  thereafter,  a 
controller  is  designed  for  the  homogeneous  half  alone  so  that  its 
output  tracks  the  difference  between  the  system  output  and  the 
output  of  the  particular  component.  The  homogeneous  half  is  a 
linear  time-invariant  system  with  desirable  convergence  prop¬ 
erties,  and  a  feed-forward  tracking  control  law  can  be  readily 
designed  for  it. 

Recall  that  the  output  of  interest  is  {t^y)dy  =  L), 

from  (7).  It  represents  the  integrated  value  of  the  wing  dihedral 
which,  in  turn,  is  a  measure  of  the  side  force  produced  by  the 
wing.  In  order  to  ensure  that  this  output  tracks  the  desired  ref¬ 
erence  signal  R{t),  the  control  signal  u{t)  [also  =  0)1  is 

designed  using  trajectory  planning  to  ensure  that  tracks 

R{t)  —  (t,  L).  In  particular,  is  approximated  via  a  polyno¬ 

mial  expansion  that  involves  powers  of  y,  and  the  coefficients 
are  solved  using  the  PDE  and  the  boundary  conditions. 

Remark  8:  The  observer-based  approach  that  is  presented 
here  can  be  used  for  a  wide  class  of  functions  F{t^y)  which, 
as  we  argued  earlier,  are  usually  nonlinear  in  ^y.  Therefore, 
it  is  difficult  to  construct  a  polynomial  expansion  for  F{t,y). 
Although  trajectory  planning  involves  a  polynomial  expansion, 
it  is  only  done  for  the  homogeneous  half,  and  consequently, 
there  is  no  need  to  obtain  a  polynomial  expansion  for  F(t,  y). 

Remark  9:  The  error  between  the  perturbation  observer  and 
the  system  is  uniformly  bounded,  and  the  bound  can  be  made 
arbitrarily  small  by  increasing  7^.  The  trajectory  planning  ap¬ 
proach  creates  a  low-pass-filtered  control  signal,  as  shown  in 
the  next  section,  and  therefore,  the  stability  and  the  robustness 
of  the  system  are  not  affected  if  a  large  value  is  chosen  for  7^ . 

Remark  10:  One  limitation  of  using  a  perturbation  observer  in 
the  present  form  is  that  it  requires  a  knowledge  of  ^(t) .  However, 
^(t)  can  be  obtained  using  an  array  of  sensors  and  spline  fitting. 
Such  an  array  of  sensors  is  light  and  cheap  enough  to  implement 
in  a  practical  setting,  and  more  so  compared  with  a  distributed 
actuation  scheme. 

D.  Trajectory  Planning  for  the  Homogeneous  Component  of 
the  Observer  (^h) 

Lor  a  reference  signal  R{t)  from  (7),  let  Rh  (t)  =  R{t)  — 
^p{t,L)  denote  the  reference  signal  that  has  to  be  tracked  by 
L),  as  shown  in  Lig.  4.  It  may  be  possible  to  use  more  than 
one  method  to  design  a  control  signal  u{t),  which  ensures  that 
the  output  of  the  homogeneous  half  tracks  the  reference  signal. 
This  freedom  results  largely  from  the  fact  that  the  homogeneous 
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component  is  an  exponentially  stable  linear  time-invariant  sys¬ 
tem  with  desirable  convergence  properties  (by  design).  In  this 
section,  we  illustrate  a  trajectory -planning  approach,  as  we  did 
for  wing  twist. 

The  term  y)  can  be  approximated  by  a  polynomial  of 
the  form 

&«.!/)  =  1;^  (59) 

i=i 

where  N  needs  to  be  chosen  to  ensure  desirable  tracking  prop¬ 
erties.  Note  that  the  coefficient  of  is  zero  because  0)  =  0. 

Substituting  (59)  into  the  dynamics  in  (52)  yields  the  set 
of  ODEs 


Vj  +  hPVj  +  =  -hVj+^  -  cib'nj+^‘  (60) 

The  boundary  conditions,  together  with  {t,  L)  =  (f),  yield 

the  following  algebraic  equations: 

Vj+2{t)U  ^  ^  ^  ^ 

j=0  j=0 

i;  aMi = (61) 

i=i 

The  value  of  N  can  be  chosen  to  ensure  good  tracking.  Note 
that  the  dynamics  are  of  fourth  order  in  y.  Therefore,  we 
need  to  choose  N  >  b,  and  the  gain  p  must  be  chosen  to  ensure 
stability  of  the  resulting  set  of  ODEs.  Eor  any  choice  of  N, 
we  have  —  4  differential  equations  (for  r]i  to  tjn-a)  and  3 
constraints,  a  total  of  —  1  equations.  The  underdetermined 
problem  may  be  resolved  by  setting  one  variable  among  r]N-3, 
r]N-2,  and  r]N-i  to  zero,  or  imposing  an  additional  constraint 
on  the  system. 

Suppose  that  we  set  A^  =  5,  i.e.,  =  ^7  =  *  •  •  =  0-  We  also 

set  7^4  (t)  =  0.^  This  yields  the  differential  equation 

fji  (t)  +  bbpyi  (t)  +  abprji  (t)  +  657)5  (t)  +  (t)  =  0.  (62) 


Constraint  (61)  can  be  solved  to  obtain  775  (t)  in  terms  of  771  (t) 
and  Rh  (t) 


vdt)  =  ^  ■  (63) 


From  (59),  it  follows  that  the  main  bending  control  input  is  given 
by 

u{t)  =  ih,v  {t,  0)  =  rii  (t)  (64) 

where  yi  (t)  is  obtained  by  integrating  the  following  ODE,  ob¬ 
tained  in  turn  by  substituting  (63)  into  (62): 


Vi  (t)  +  h 


f]i  (t)  +  ab 


{bbRh{t)  +  abRh{t)). 


(65) 


^Motivation:  if  is  a  constant,  then  774  =  0  is  in  the  steady-state  solution. 
Therefore,  this  approximation  works,  at  least,  for  a  class  of  slow  time- varying 
signals. 


Note  that  the  gain  p  has  to  be  chosen  so  that  p  >  12/ for 
775  (t)  to  be  stable.  This  completes  the  design  of  a  trajectory 
planning-based  tracking  controller  for  bending. 


VI.  Simulations 


Simulations  are  carried  out  in  MATLAB  by  using  a  Galerkin- 
based  approach  to  convert  the  PDE  system  into  ODEs.  The 
Galerkin  truncation  is  not  used  as  a  basis  for  control  law  design; 
therefore,  no  danger  of  “spillover  instability”  arises.  The  twist 
6>(t,  y)  is  expressed  as  a  weighted  sum  of  basis  functions  (pi  (7/), 
7  =  1,  2, . . . ,  77  and  ^p{y) 

n 

0{t,  y)  =  s{t)'tp{y)  +  '^Vi  (y) 

i=l 

cPi{0)=4{L)=0  (66) 

where  s{t)  is  the  boundary  control  input.  If  boundary  control  is 
applied  at  the  wing  root,  then  7/)  (7/)  has  to  be  chosen  to  satisfy 
'ip'{L)  =  0  and  7/)(0)  =  1.  On  the  other  hand,  if  the  boundary 
control  is  applied  at  the  wing  tip,  we  choose  to  satisfy  7/)(0)  = 
0  and  7/)'(I/)  =  1. 

The  PDE  in  (32)  can  be  rewritten  as 

n  n 

"k  ^  ^^yi (7/)  ~ 

i=l  i=l 

-  4>"{y){bs  +  as) 

=  M  i^{y)s  +  y  i()4>i  (y)  j  •  (67) 

Using  Galerkin’s  method,  (67)  is  converted  to  a  set  of  ODEs 


c{s{t)  —  Ms{t))  +  d{bs  -h  as{t)) 

+  -  Mri{t))  +  [B]{bfi{t)  +  ari{t))  =  0  (68) 

where 

C=  f  'ip{y)(l){y)dy,  d  =  - 
Jo 

[^]  =  Y  cjxfdy,  [B]  =  -  Y  cf>{cp")^dy.  (69) 
Jo  Jo 

The  reader  will  recall  that  t/)  is  a  scalar,  and  (p  =  [0i ,  02 ,  •  •  •  ^ 

is  a  vector.  The  control  s{t)  is  expressed  similarly  in  terms  of 
0,  0,  and  y  to  obtain  a  set  of  ODEs,  which  are  simulated  to 
approximate  the  response  of  the  system. 


/ 


i’"{y)Hy)dy 


A.  Root  Control  of  Wing  Twist 

Eig.  5  demonstrates  the  regulation  of  twist  dynamics  us¬ 
ing  the  backstepping  controller  that  is  derived  in  (31),  with 
the  transformation  in  (12)  and  (17).  The  value  of  M/a  was 
set  to  8,  where  a  =  GJ /Ip.  A  value  of  p  =  4  yielded  an 
unstable  response,  while  the  response  was  stable  for  p  =  8. 
Recall  the  following  condition  for  stability  with  L  =  1:  p  > 
M/a  —  7r^/4?^5. 5.  The  backstepping  controller  works  even 
when  M{y)  =  M(1  —  p^)  is  used  (to  mimic  an  elliptical  lift 
distribution  over  the  wing)  instead  of  a  constant  M{y)  =  M. 
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Fig.  5.  Regulation  of  the  twist  dynamics  using  the  backstepping  controller  in 
(31).  The  plots  were  obtained  for  M/a  =  8,  while  p  was  increased  to  ensure 
stability.  Each  plot  is  a  collection  of  appropriately  chosen  snapshots  in  the 
interval  0  <  t  <  30  s,  with  the  lines  getting  darker  with  time.  In  plots  (b)  and 
(c),  6{t,  y)  converges  to  zero  exponentially,  (a)  Unstable  response  with  p  =  4. 
(b)  Stable  response  with  p  =  8.  (c)  Stable  response  with  p  =  8  for  spatially 
varying  M(y). 
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The  backstepping  controller  can  be  added  on  top  of  a  tracking 
controller. 

Each  plot  in  Fig.  5  is  a  collection  of  snapshots  that  are  ap¬ 
propriately  chosen  in  the  interval  0  <  t  <  30  s,  with  the  lines 
getting  darker  with  time.  The  system  state  0{t^y)  converges  to 
zero  exponentially  when  p  =  8. 


Fig.  6.  Twist  profile  of  the  wing  as  a  function  of  time  when  the  the  adaptive 
controller  in  (41)  is  applied  at  the  wing  tip.  Three  cases  have  been  examined 
pL 

here,  with  6>(t,  =  0.05  as  the  desired  output.  The  first  three  plots  are 

appropriately  chosen  snapshots  in  the  interval  0  <  t  <  30  s,  with  lines  get¬ 
ting  darker  with  time.  The  error  metric  e{t)  converges  to  zero  exponentially, 
(a)  M  unknown,  constant,  open  loop  stable,  (b)  Unknown  M,  constant, 
open  loop  unstable,  (c)  Spatially  varying,  unknown  M  (y),  unstable  dynamics, 
(d)  Time  history  of  e(t) 


B.  Tip  Control  of  Wing  Twist 

Fig.  6  shows  simulation  results  for  the  twist  dynamics,  when 
the  actuator  is  located  at  the  wing  tip.  The  value  of  M/a  was 
chosen  so  that  stability  is  assured  without  the  need  for  a  ded¬ 
icated  stabilizing  controller.  The  first  plot  was  obtained  for  a 
system,  where  the  aerodynamics  were  assumed  to  be  linear  but 
unknown.  The  second  plot  was  assumed  linear  unknown  aero¬ 
dynamics  such  that  the  system  was  unstable  in  the  open  loop, 
but  within  the  enhanced  stability  margin  that  is  described  in 
Section  IV-B.  The  third  plot  (bottom  row)  was  obtained  for  the 
case  where  the  aerodynamics  were  additionally  spatially  vary¬ 


ing.  The  fourth  plot  depicts  the  time  histories  of  the  tracking 
error  e(t)  and  the  control  signal  u{t),  respectively.  In  all  cases, 
the  twist  amplitude  converges  to  the  steady-state  value  with 
satisfactory  transients.  The  error  metric  (5)  converges  to  zero 
exponentially. 

C.  Comparison  of  Finite-State  Approximation-Based  and 
Trajectory  Planning-Based  Controllers  for  Bending 

Using  Galerkin’s  method,  the  beam  bending  dynamics  (46), 
with  the  control  input  u{f)  and  its  derivative  u{f)  as  additional 
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Fig.  7.  Comparison  of  controllers  derived  from  finite-state  approximations  of 
the  open  loop  and  the  homogeneous  component,  respectively,  with  a  controller 
obtained  using  trajectory  planning  [41].  (a)  Traditional  ODE  control,  (b)  PDE 
perturbation  observer  with  trajectory  planning. 

states,  can  be  recast  into  the  form 

X  =  Ax  +  Bu  +  f{x),y  =  Cx  (70) 

where  u{t)  equals  the  highest  derivative  of  u{t),  i.e.,  u{t)  = 
u{t).  The  system  can  be  stabilized  if  (A,  B)  is  stabilizable.  In 
particular,  vectors  B  and  C  depend  strongly  on  the  choice  of  the 
basis  functions.  The  unmatched  nonlinearity  f{x)  can  affect  the 
tracking  error  as  well  as  the  stability  of  the  closed-loop  system 
if  not  compensated  adequately.  The  unmatched  nonlinearity  can 
be  tackled  using  backstepping,  or  a  more  direct  approach  such 
as  in  [23,  Ch.  3,pp.  140-174]. 

Fig.  7  compares  the  performance  of  two  controllers.  The  first 
is  based  on  [23]  and  designed  using  a  finite-state  approximation 
of  the  beam,  and  the  other  is  based  on  the  perturbation  observer 
and  trajectory  planning  approach  in  Section  V.  The  beam  is 
assumed  to  be  loaded  uniformly,  but  for  the  purpose  of  control 
design,  this  load  is  treated  as  an  unknown.  This  simple  loading 
condition  is  enough  to  create  an  unmatched  uncertainty  in  the 
finite- state  approximation. 

Fig.  7(a)  shows  the  performance  of  a  controller,  which  is  de¬ 
signed  using  the  finite-state  approximation.  The  tracking  error 
is  less  than  5%  of  the  signal  amplitude,  and  there  is  a  phase  lag 
in  the  system  output  when  compared  with  the  reference  signal. 
The  tracking  error  as  well  as  the  phase  lag  arises  due  to  the  use 
of  a  low-pass  filter  in  the  generation  of  the  control  signal.  On  the 
other  hand.  Fig.  7(b)  [41]  shows  the  tracking  performance  of  the 
proposed  controller  which  is  derived  from  trajectory  planning 
using  a  seventh  degree  polynomial  and  the  perturbation  observer 
in  Section  V-C.  The  amplitude  of  the  tracking  error  is  consid¬ 
erably  smaller  than  the  previous  approach,  and  the  phase  lag  is 
almost  eliminated.  For  the  system  that  is  considered  in  this  il¬ 
lustration,  the  proposed  controller  outperforms  the  conventional 
ODE-based  controller.  Note  that  the  trajectory  planning-based 


Fig.  8.  Experimental  setup  showing  the  beam  and  power  supply. 

control  signal,  which  is  generated  using  (62)  and  (64),  is  it¬ 
self  a  low-pass  filter,  albeit  one  that  ensures,  by  design,  a  small 
tracking  error. 

VII.  Experimental  Results 

In  this  section,  we  describe  some  successful  experiments, 
which  were  performed  to  validate  the  efficacy  of  the  pertur¬ 
bation  observer-based  controller  in  Section  V.  Fig.  8  shows 
the  experimental  setup.  A  long  thin  beam  with  a  rectangular 
cross  section  is  utilized  for  experiments.  It  is  rigidly  clamped 
at  one  end,  while  the  other  end  is  free  to  be  controlled.  Two 
small  permanent  magnets  attached  to  the  tip  are  acted  on  by  an 
electromagnet  whose  magnetic  field  can  be  controlled  by  spec¬ 
ifying  a  current  and  varying  the  voltage  using  a  programmable 
power  supply  (Agilent  Technologies  E3642A)  that  interfaces 
with  MATLAB.  The  interaction  between  the  permanent  mag¬ 
nets  and  the  magnetic  field  of  the  electromagnet  produces  a 
force  that  acts  on  the  beam  tip. 

Retroreflective  markers  are  placed  at  seven  specific  points 
along  the  length  of  the  beam,  and  the  VICON  motion  capture 
system  is  used  to  track  the  coordinates  of  those  markers.  Bending 
displacement  is  determined  by  calculating  the  distance  between 
the  markers  on  the  beam  in  the  deformed  state,  and  the  coor¬ 
dinates  of  the  corresponding  markers  in  the  undeformed  state. 
The  bending  displacement  is  interpolated  from  the  data  that  are 
obtained  from  the  VICON  system  to  obtain  the  deformation  pro¬ 
file.  This  information  is  used  to  calculate  the  tip  force,  which  is 
the  control  variable. 

The  electromagnet  is  characterized  to  determine  a  relation¬ 
ship  between  the  required  voltage  and  the  commanded  tip  force. 
The  tip  deflection  ^tip  of  a  cantilever  beam  as  a  function  of  the 
applied  tip  force,  (=  El^^yyy  (t,  L)),  is  given  by 

Ctip=C(i,U  =  f||^  ^  u{t)=^yyy{t,L)  =  ^.  (71) 

The  electromagnet  is  calibrated  by  measuring  the  steady-state 
tip  deflection  as  a  function  of  the  input  voltage  so  that 

^tip  «  1.33  X  (72) 

where  W  denotes  the  voltage  of  the  signal,  which  is  sent  to  the 
electromagnet.  Combining  (71)  and  (72)  with  the  fact  that  the 
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Fig.  9.  Time  histories  of  the  tip  displacement  (^(t,  L)),  when  a  sinusoidal 
reference  signal  was  sent  to  the  system.  An  impulse-like  disturbance  was  also 
administered  to  the  beam.  The  controller  was  disabled  in  the  first  case  to  provide 
a  comparison  between  the  open-  and  closed-loop  systems,  (a)  Open  loop  with 
an  impulse  perturbation,  (b)  PDF  boundary  control  with  impulse  perturbation. 

control  input  u{t)  =  ^yyy{t,  L)  =  Ftip/E'/^  results  in 
4x1 

u(t)  =  — ®  ^  W  =  39.67L(u(t)f/\  (73) 

The  control  law  u(t)  is  obtained  from  an  equation,  which  is 
analogous  to  (64)  [19] 

^8Rh(t)  ^  £  u  T\ 

u{t)  —  ^  ,  Rfiit)  —  R\t) 


VIIL  Conclusion 

This  paper  has  introduced  a  boundary  control  formulation  for 
a  class  of  PDEs  whose  output  consists  of  a  spatial  integral  of 
the  state  variable.  Although  the  primary  focus  of  this  paper  was 
on  the  control  of  flexible  wings,  this  boundary  control  problem 
is  directly  relevant  to  the  control  of  flexible  continuum  robotic 
arms  described  by  PDEs. 

In  order  to  control  wing  twist  using  a  root-based  actuator, 
a  tracking  controller,  which  was  based  on  PDE  backstepping 
was  designed  to  handle  the  output  that  has  an  infinite  relative 
degree  with  respect  to  the  root  actuation.  On  the  other  hand, 
when  twisting  moment  at  the  wing  tip  was  the  control  input, 
it  was  shown  that  the  system  had  a  relative  degree  of  2  with 
respect  to  the  input-output  combination,  which  facilitated  the 
design  of  a  robust  adaptive  tracking  controller.  It  was  shown  that 
the  tracking  controller  alone  could  yield  a  16-fold  improvement 
in  the  stability  margin  of  the  system. 

The  force  distribution  that  drives  bending  is  far  too  com¬ 
plicated  to  be  handled  by  backstepping  alone,  which  motivated 
the  development  of  a  new  perturbation-observer-based  approach 
that  is  presented  in  Section  V.  The  two-stage-observer-based 
tracking  controller  was  designed  to  ensure  that  the  wing  tip 
tracked  the  desired  reference  profile  with  the  bending  PDE.  The 
first  component  of  the  observer  compensated  for  the  nonlinear¬ 
ities  and  unmodeled  dynamics,  while  the  other  was  used  exclu¬ 
sively  for  deriving  the  control  signal.  The  proposed  perturbation- 
observer-based  controller  was  implemented  in  an  experimental 
setup,  which  allowed  us  to  control  the  bending  displacement 
of  the  tip  of  a  beam  by  applying  a  force  at  the  wing  tip.  The 
controller  performed  as  expected,  and  demonstrated  reference 
signal  tracking  as  well  as  disturbance  rejection. 
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The  gain  was  set  to  p  =  500  for  the  experiments.  Pig.  9  shows 
the  time  history  of  the  tip  displacement  for  two  cases.  In  both 
cases,  the  system  was  given  an  impulsive  disturbance  midway 
through  the  experiments  at  t  =  17  s,  in  addition  to  persistent 
external  disturbances  in  form  of  a  mild  breeze  in  the  experimen¬ 
tal  area.  The  open-loop  response  is  noisy  due  to  the  persistent 
external  disturbances.  The  amplitude  takes  a  very  long  time  to 
converge  after  the  impulsive  perturbation  at  17  s.  In  comparison, 
the  closed-loop  response  shows  a  negligible  tracking  error  and 
rapid  convergence  to  the  reference  signal  after  the  impulsive 
perturbation.  Additional  experimental  results,  including  those 
for  root  control  of  beam  bending,  can  be  found  in  [19]  and  [42]. 
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Dyadic  Perturbation  Observer  Framework  for 
Control  of  a  Class  of  Nonlinear  PDE/ODE  Systems 

Aditya  A.  Paranjape  and  Soon-Jo  Chung 


Abstract — This  paper  presents  the  general  theory  of  the  dyadic 
perturbation  observer  (DPO)  framework  as  a  generic  method  for 
controlling  a  class  of  systems  described  by  partial  and/or  ordinary 
differential  equations.  The  method  is  particularly  applicable  to 
boundary  control  problems  for  systems  described  by  partial 
differential  equations  (PDEs).  Conditions  for  closed-loop  stability 
and  robustness  are  derived  using  the  small  gain  theorem,  and  the 
results  are  further  specialised  for  finite  dimensional  systems. 


1.  Introduction 

The  Dyadic  Perturbation  Observer  (DPO)  framework  was 
first  presented,  under  the  title  of  perturbation  observer-based 
control,  in  [12]  although  its  closed-loop  stability  was  not 
explicitly  proved.  The  method  was  demonstrated  successfully 
in  experiments  involving  control  of  beam  bending  in  [11], 
[12].  The  objective  of  this  note  is  to  present  the  general 
theory  of  the  DPO  framework  as  a  generic  method  for  a 
class  of  systems  described  by  ordinary  and  partial  differential 
equations  (ODEs  and  PDEs),  including  combinations  of  both, 
and  derive  conditions  for  closed-loop  stability  and  robustness. 

Systems  described  by  a  combination  of  ODEs  and  PDEs 
arise  naturally  in  engineering  applications,  such  as  fiexible 
aircraft  wings  (fiexible  wing  structure  combined  with  the  rigid 
body  aircraft  dynamics),  robotic  surgical  systems  (a  multi- 
segmented  fiexible  robotic  arm),  temperature  control  systems 
(heat  diffusion  and  mass  fiow  of  air),  etc.  There  are  several 
well-  established  families  of  methods  for  designing  controllers 
for  finite  dimensional  systems,  such  as  dynamic  inversion, 
backstepping,  gain  scheduling,  and  Lyapunov  function-based 
approaches.  The  important  feature  of  all  these  methods  is  that 
they  apply  to  any  given  class  of  ODEs,  although  it  may  not 
always  be  convenient  to  apply  one  particular  method  to  a  given 
system.  In  comparison,  systems  described  by  nonlinear  PDEs 
are  controlled  either  through  their  ODE  approximations  [1], 
[2]  or  using  methods  which  leave  the  PDE  intact  and  yield 
closed-loop  stability  and  performance  guarantees,  but  tend  to 
require  information  about  the  exact  structure  of  the  PDE  [6], 
[8],  [9],  [13]. 

The  DPO  framework  is  designed  primarily  for  the  following 
class  of  PDE  systems,  although  it  applies  equally  well  to  ODE 
systems: 

Wt{t,  x)  =  =  Aw{t,  x)  +  f{t,  X,  w{t,  x))  (1) 
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Champaign.  Email:  {paranja2,  sjchung}@ illinois.edu.  This  research  was 
supported  by  NSF  (IIS-1253758)  and  ARO  (W91  INF- 10- 1-0296). 


where  w{t,x)  is  the  state  of  the  system,  x  is  the  spatial 
coordinate,  and  w  denotes  the  vector  of  w  and  its  partial  spatial 
derivatives.  The  problem  formulation  will  be  made  precise 
later  in  the  paper.  Note  that  the  operator  ^  is  a  linear,  time- 
invariant  operator.  The  nonlinearities  are  all  captured  by  the 
function  /(•).  We  denote  by  Wx  the  partial  derivative  dwjdx, 
additional  subscripts  denote  higher  order  partial  derivatives 
with  respect  to  those  variables,  and  Wx{t,x)  is  defined  like¬ 
wise. 

Example  1.  In  the  most  simplified  form,  the  bending  dynam¬ 
ics  of  a  wing  attached  to  a  moving  aircraft  are  given  by  [10] 

El  {v^txxxx  E  ^  xxxx  )=F{x,i,a,V),  (2) 

0)  =  0)  =  L)  =  ^xxxit,  L)  =  0 

q=  Y  J  +  Mtaii 

where  >  0  is  the  coefficient  of  Kelvin- Voigt  damping, 
while  denotes  the  nonlinear  aerodynamic  force 

terms  on  the  wing  per  unit  span.  The  dynamics  of  the  angle 
of  attack  a  and  the  pitch  rate  q  capture  the  rigid  body 
motion  of  the  aircraft  at  speed  V.  The  term  Mtaii  denotes 
the  control  moment  from  the  horizontal  tail  of  the  aircraft, 
Xa  is  a  reference  length  which  measures  the  distance  between 
the  wing  and  the  center  of  mass  of  the  aircraft,  and  Iq  isthe 
principal  moment  of  inertia  about  the  pitching  axis. 

The  left-hand  side  of  (2)  captures  the  elastic  behavior  of 
the  wing.  Moreover,  for  a  wing  immersed  in  an  external  fiow 
field,  the  stack  ^  is  given  by  ^  =  [^,  ^x]^  although  other 

spatio-temporal  derivatives  may  be  included  for  completeness. 
Control  design  for  this  problem,  with  treated  as  the 

control  input,  has  been  addressed  in  [12],  while  the  stability 
of  the  closed  loop  is  addressed  in  this  paper. 

The  DPO  framework  decomposes  the  system  in  (1)  into 
two  halves  for  the  purpose  of  control  design:  one  half  (called 
the  particular  half)  only  accommodates  the  nonlinearity  (/(•) 
in  (1)),  while  the  other  half  (called  the  homogeneous  half) 
accommodates  only  the  control  input  (which  could  be  a 
boundary  condition  of  a  PDE).  The  DPO  architecture  has  been 
illustrated  in  Pig.  1.  The  control  signal  is  designed  to  ensure 
that  the  output  of  the  homogeneous  half  tracks  the  desired 
reference  signal  minus  the  output  of  the  particular  half,  thereby 
ensuring  that  the  output  of  the  two  halves  put  together  tracks 
the  reference  signal.  The  stability  of  the  closed  loop  is  verified 
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explicitly  using  the  small  gain  theorem.  Note  that  the  observer 
is  used  (in  line  with  the  certainty  equivalence  principle)  to 
identify  the  particular  and  homogeneous  state  variables.  By 
nature  of  its  design,  the  DPO  framework  is  inherently  adaptive 
in  nature  so  that  it  can  accommodate  modeling  and  parametric 
uncertainties,  as  well  as  external  disturbances  (e.g.,  see  [12] 
for  the  use  of  projection-based  adaptive  control  along  with  the 
DPO). 


Fig.  1.  A  block  diagram  of  the  DPO  framework,  with  the  subscripts  p  and 
h  denoting  the  particular  and  homogeneous  components,  respectively. 

The  salient  features  of  the  DPO  method  are: 

1)  The  DPO  control  design  is  dimension-independent,  i.e., 
it  can  be  applied  to  finite  as  well  as  infinite  dimensional 
systems  (or  a  combination  thereof). 

2)  The  DPO  does  not  prescribe  any  particular  control 
design  for  the  homogeneous  half;  any  control  signal  that 
guarantees  tracking  is  admissible,  including  nonlinear 
control  signals,  as  long  as  the  resulting  closed-loop 
system  satisfies  a  small-gain  condition. 

The  paper  is  organized  as  follows.  Mathematical  prelim¬ 
inaries  are  recapitulated  in  Sec.  II,  followed  by  the  problem 
formulation  in  Sec.  III.  The  DPO  control  architecture  has  been 
presented  in  Sec.  IV,  while  a  finite  dimensional  analogue  is 
discussed  briefiy  in  Sec.  V.  Simulation  results  illustrating  the 
DPO  framework  are  presented  in  Sec.  VI. 

II.  Preliminaries 

Let  Z  be  a  Hilbert  space  with  the  usual  inner  product 
defined  by  (zi,  Z2)  =  /q  zjz2dy  for  any  2^1,  ^2  G  Z,  where 
L  denotes  the  total  span  of  the  spatial  dimension.  We  are 
concerned  with  variables  w{t,x)  G  M  x  Z,  with  x  G  [0,  L]. 
We  define  the  space  W  =  IR+  x  Z,  so  that  w{t,  x)  G  W.  The 
space  W  is  a  Banach  space  with  the  norm  defined  as  follows. 


Definition  1.  Given  w(t^x)  G  W,  we  define  the  following 
norms: 


Definition  2.  Given  an  operator  ^  :  W  i-A  W  denote  its 
induced  norm  by  \\A\\i. 

For  a  time-varying  signal  p{t)  G  M,  we  define  its  jCi 
norm  as  \\p\\ci  =  It  can  be  checked  that 

ll^lli  =  ||crmax,x(^(t))IUi,  where  crmax,x(^(^))  IS  a  time- 
varying  signal  and  the  maximum  singular  value  is  taken  over 
all  values  of  x  at  each  instant  in  time.  It  must  be  noted  that 
the  jCi  norm  is  taken  with  respect  to  time. 

Definition  3.  Given  w{t^x)  G  W,  let  w(t,x)  denote  a  vector 
of  w  and  its  partial  spatial  derivatives  with  orders  less  than 
that  of  the  highest  order  in  the  operator  A.  This  is  a  fairly 
straight-forward  extension  and  allows  for  accommodating 
forcing  functions  that  depend  on  the  partial  derivatives  of  w. 

Assumption  1.  We  set  all  initial  conditions  to  zero.  This  does 
not  alter  the  fundamental  nature  of  the  result  because  it  only 
introduces  exponentially  decaying  terms. 

Next,  we  state  a  variation-of-constants-like  formula  for 
infinite  dimensional  systems  (Lemma  3.1.2,  [3]). 

Lemma  1.  Consider  a  system  w  =  Aw  +  /(t,  x,  w),  w{t)  = 
Wo  G  D{A)  C  W,  where  A  is  the  infinitesimal  generator  of  a 
semigroup  T{t)  and  D{A)  denotes  the  domain  of  A.  Suppose 
that  /(•)  is  a  smooth  function  of  its  arguments.  Then,  w{t)  is 
given  by 

w{t)  =T{t)wo^  [  T{t-T)f{r,x,w{r))dr 
Jo 

where  wo  =  w{0,  •),  which  is  set  to  zero  by  Assumption  1. 

Definition  4.  Given  a  semi-group  T (t),  we  define  the  operator 
T (t)x  :W  so  that  V  g(t,  x)  G  W  and  V  t  >  0, 

T{f)'kq{t,x)  =  f  T{t  -  r)q{r,x)  dr 
Jo 

We  further  define  T'k  =  limt^oo  (^nd  the  induced  norm 

lir*||i  =  sup||r(t)*||i 

t>0 

In  order  to  help  establish  point- wise  boundedness  in  space 
(as  against  integral  boundedness  in  the  sense  of  Z),  we 
introduce  Agmon’s  inequality  (Lemma  2.4,  [8]). 

Lemma  2  (Agmon’s  inequality).  Given  w{t,x)  satisfying 
||u;(t)||^  <  00  and  ||u;a^(t)||z  <  00,  we  have 

max  \w{t,x)\‘^  <  \w{t,0)\‘^  +  2\\w{t)\\z\\w^{t)\\z 

0<x<L 

Agmon’s  inequality  allows  us  to  bound  the  £00  norm  of 
signals  in  the  spatial  domain  provided  the  spatial  Z  norms  of 
the  signals  and  their  first  derivatives  are  bounded. 

Definition  5.  Given  p{t)  G  with  components  Pi{t)  (1  < 
i  <  n),  we  define 

IbWIloo  =  \pi{t)\ 

l<i<n 

Ibikoo  =  esssup  ||p(t)||oo,  ||p||£^,x  =  ess  sup  ||p(t)||oo 

t>0  0<t<r 

^  denote  p  G  £^. 
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Definition  6.  The  Ci  norm  of  a  linear  operator  T  : 
is  defined  as  ||7'||£i  =  sup||p||^^=i  II-^pIUoo^  P  e 


III.  Problem  Formulation 

We  consider  a  class  of  systems  described  by  PDFs  of  the 
form 


Wt  =  +  X,  w),  Bw  =  u{t),  Nw  =  0,  y{t)  =  Cw  (4) 

where  w  =  w{t,x);  u{t),  y{t)  eR,  and  A  is  the  infinitesimal 
generator  of  a  continuous  semigroup  T{t).  The  right-hand 
side,  /(t,x,  w),  is  a  nonlinear  function  of  its  arguments.  The 
operators  ff  and  B  capture,  respectively,  the  homogeneous 
boundary  conditions  and  the  boundary  control  input.  The 
control  objective  is  to  ensure  that  the  toutput  y{t)  tracks  a 
reference  signal  r{t)  G  M. 

Assumption  2.  The  semigroup  T{t)  (whose  infinitesimal 
generator  is  A  in  (4))  satisfies  ||T(t)||i  <  >  0, 

where  uj  >{)  is  a  constant.  Moreover,  \\T ^  ||i  is  bounded. 

Thus,  the  system  dynamics  can  be  viewed  as  the  sum 
of  a  linear,  exponentially  stable,  well-posed  operator  and  an 
external  nonlinear  forcing  term. 

Assumption  3.  For  every  p  >  0,  there  exist  positive  constants 
ui(p)  and  U2{p)  such  that  if  ||u;||w,r  <  pfar  some  r  >  0,  then 

||/(t,X,  W)||w,r  <  Z/l(p)||u;||w,r  +  Z^2(p), 

In  general,  ui(p)  and  U2{p)  are  class  IC  functions  of  p. 


The  above  assumption  essentially  implies  that  the  spatial 
derivatives  of  w{t,x)  enter  through  bounded  functions.  As 
an  illustration,  in  Example  1,  the  spatial  derivative 
enters  through  the  trigonometric  functions  cos(^x(t,  x))  and 
sm(^^(t,x)). 


Assumption  4.  The  output  operator  C  in  (4)  is  bounded;  i.e., 
WvWjCoo  =  ll^'^(^5^)IUoo  ^  ^ll'^llw  for  some  constant  K. 


Example  2.  Let  y{t)  =  w{t,x)dx.  Then,  |^(t)|  = 

f^w{t,x)dy  <  VL\\w{t)\\z  IIpIUoo  ^  \/i||w||w. 

Note  that  the  inequality  is  a  consequence  of  w{t,x)  lying  in 
a  Hilbert  space.  Assumption  4  is  relaxed  in  Sec.  IV-C.  An 
example  of  an  unbounded  operator  is  the  value  of  w{t,x)  at 
a  boundary,  e.g.,  Cw{t,x)  = 


IV.  Control  Design  using  DPO 
A.  Design  of  the  Perturbation  Observers 

We  use  the  operator  "  to  denote  observer  states,  and  the 
subscripts  p  and  h  to  denote  states  of  the  particular  and  the 
homogeneous  halves,  respectively.  The  dynamics  of  the  two 
halves  are  given  by 

=  Awp  +  f{t,  X,  w),  Bwp  =  Afwh  =  0,  yp  =  Cwp  (5) 

Wh,t  =  Awh,  Bwh  =  u{f),  ffwh  =  0,  yh=  Cwh  (6) 

Recall  that  the  subscript  t  denotes  the  partial  derivative  d/dt. 
Using  the  projection  operator  to  design  /(•),  such  as  in  [12], 
we  can  show  the  following  result. 


Lemma  3.  ^  ||u;||w  <  p,  then  there  exist  constants  no  = 
nQ{p)  and  =  ni(p)  such  that  ||u)p||w  <  /^o||'^||w  + 

Proof.  The  proof  follows  by  applying  Lemma  1  to  (5)  and 
computing  the  norm  of  Wp.  ■ 

We  define  the  observation  error  w  =  Wp  Wh  —  yo.  Then, 
from  (4),  (5),  and  (6),  we  conclude  that  the  dynamics  of  w 
are  given  by 

wt=  AwF  f{t,  X,  w)  -  f{t,  X,  w), 

Bw  =  Mw  =  0,  y{f)  =  x)  (7) 

Lemma  4.  It  is  possible  to  ensure  that  is  uniformly 

bounded  V  t  >  0,  and  the  bound  can  be  made  arbitrarily  small. 
It  follows  as  a  corollary  that  y{t)  is  uniformly  bounded  and 
can  be  made  arbitrarily  small. 

Proof.  The  boundedness  of  \\w\\z  is  proved  in  [12]  (Theo¬ 
rem  4),  while  that  of  y{t)  follows  from  Assumption  4.  ■ 


B.  DPO  Control  Synthesis  and  Closed-Loop  Stability  Analysis 

We  design  the  control  signal  u{t)  to  ensure  that  the  output 
of  the  homogeneous  half,  yh{t)  in  (6),  tracks  a  reference  signal 
r(t)  —yp{t),  where  r{t)  is  the  reference  signal  for  the  original 
system  (4),  where  yp  is  the  output  of  the  particular  half  (5). 
Lemma  4  would  then  ensure  that  the  output  y{t)  of  the  system 
(4)  tracks  the  reference  signal  r(t),  as  desired.  We  can  write 
the  input-output  dynamics  of  the  linear,  exponentially  stable 
homogeneous  half  in  (6)  in  the  Laplace  domain  [4]:  yh{s)  = 
Gc{s)u{s),  where  the  transfer  function  Gc{s)  depends  on  A 
and  the  boundary  conditions,  as  explained  in  [4].  We  design 
u{t)  in  the  Laplace  domain  as  follows: 

u{s)  =  H{s){r{s)  -  yp{s))  =  H{s){r{s)  -  Cwp{s))  (8) 

The  ideal  choice  for  77(5)  would  then  be  1/Gc{s)  which  is,  in 
general,  neither  proper  nor  stable.  Therefore,  we  choose  H (s) 
to  satisfy  77(0)Gc(0)  =  1  and  impose  further  conditions  for 
closed-loop  stability  (see  Theorem  1).  Interestingly,  a  poly¬ 
nomial  approximation  of  the  state  Wh  yields  a  low-pass  filter 
77(5)  which  ensures  that  u{t)  is  sufficiently  differentiable,  and 
all  derivatives  are  bounded  [11],  [12]. 

The  next  result  asserts  the  boundedness  of  control  inputs. 

Lemma  5.  ^  ||u;||w  <  p,  then  the  control  law  in  (8) 
ensures  that  there  exist  constants  Sir  =  ^ir(p^  ^{^))’  ^iw  = 
Siw{p,  H{s)),  and  Siu  =  Siu{p,  77(5))  such  that  < 

(^irlklUco  +  +  5iu,  where  and 

=  u{t). 

Proof:  Lrom  (8),  we  get  ||^||£^  <  ||77(5)||£, (||r||/:^  + 
\\Cwp\\jc^).  The  proof  follows  from  Assumption  4  and 
Lemma  3.  ■ 

We  will  prove  the  bounded-input-bounded-output  (BIBO) 
stability  of  the  closed-loop  (in  the  sense  of  £00  boundedness 
of  signals)  using  the  small  gain  theorem.  The  first  result  asserts 
that  the  boundary  control  system  in  (4)  can  be  recast  into  a 
PDE  system  with  homogeneous  boundary  conditions  (Theorem 
3.3.3,  [3]). 
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Lemma  6.  The  coordinate  transform  v  =  w  —  jSu,  where  the 
operator  jS  satisfies  Bfdu  =  u  and  M (du  =  0,  transforms  the 
system  (4)  into  the  form 

V  =  Av  +  A(du  —  l3u-\-  f(t,  y,  w),  Bv  =  J\fv  =  0  (9) 

Consider  the  system  in  (9),  and  suppose  ||u;(t)||z  =  p  for 
some  time  t  >  0  and  furthermore,  ||u;(t)||z  <  p  Vr  <  t. 
Using  Lemma  1,  Assumption  3  and  Lemma  5,  we  get 

||'^||w,t  <  \\T ^  ||i(^o||'^||w,t  +  ^ilklUoo  +  ^2) 

where  ^0  =  z^o(p)  +  \\Af3\\i5ow  +  =  ||^/3||i^0r  + 

||/3||i^i^,  and  62  =  z/i(p)  +  WA^WAou  +  W^WAiu-  Since  w  = 
V  +  (5u,  using  the  triangle  inequality  and  Lemma  5,  we  get 


||'^||w,t  <  \\T  'k  ||i(^o||'^||w,t  +  ^llklUoo  +  ^2) 
+  ||/5||i(^0w  ||'^||w,t  +  ^OrlklUoo  +  ^Ou) 
Ai||r||£^  +  A2 


1- Ao 


(10) 


^0  —  \\T  k:  ll^^o  +  ||/3||i^Ow5  —  ||T  'k  ||^^i  +  ||/3||i^Or5 

A2  =  ||r^||i^2  +  ||/3||i^0n  (11) 


We  prove  the  following  result. 


Theorem  1.  The  closed-loop  system  (4),  (5),  (6),  and  (8)  is 
BIBO  stable  in  the  sense  of  Coo  If  there  exists  a  stable,  strictly 
proper  H{s)  and  a  constant  p  such  that 


Ai||r||/:^  +  A2 

1  —  Ao 


<  p-e 


(12) 


where  Ai's  are  defined  in  (11),  and  e  >  0  A  arbitrarily  small 


Proof.  Suppose  that  ||u;||w,t  =  ||^(f)||z  =  P  for  some  t  >  0. 
Then,  from  (10)  and  (12),  it  follows  that  ||u;||w,t  ^  p  —  ^  <  P^ 
which  contradicts  the  assumption  that  ||u;||w,t  =  p,  and  shows 
that  ||u;||w  <  P  Vt.  This  completes  the  proof.  ■ 


C.  Redressing  Assumption  4  and  Generalization 

It  is  possible  to  relax  Assumption  4  by  deriving  a  stability 
condition  of  the  form  in  Theorem  1  for  a  system  consisting 
of  the  dynamics  of  w  rather  than  just  w.  Recall  that  w  is  the 
vector  of  w  and  its  partial  spatial  derivatives.  The  resulting 
stability  condition  is  cumbersome  and  difficult  to  verify  in  a 
practical  setting;  its  representation  in  a  more  tractable  form  is 
an  open  problem.  It  must  be  noted  that  the  w  dynamics  are 
constructed  purely  as  an  analytical  tool  to  verify  the  stability 
of  the  closed  loop  system  (4),  (5),  (6),  and  (8). 

We  start  by  showing  a  number  of  properties.  Given  the 
system  (4),  we  can  differentiate  the  state  variable  w{t,  x)  with 
respect  to  x  and  recast  it  into  the  form 

w  =  ^(w)  W  +  /  (t,  X,  w) 

Sw  =  u{t),  Jfw  =  0,  y{t)  =  Cw  (13) 

Example  3.  Consider  the  PDE  system  +  ^xx  =  sin(^a^)  + 
x^,  with  the  boundary  control  ^(t,0)  =  u{t).  The  array  ^  = 
[^,  From  the  above  equation,  we  have  that 

^xt  3“  ^XXX  —  cos{Cx)C  XX  2  X, 


SO  that,  using  the  notation  from  (13), 

^(0  = 


■9*2 


0  cos{^x) 
f{t,x,0  =  [sin(^a:) 


0 

_2_  _  _2!_ 
dx  dx‘^ 


Note  that  the  operator  A{^)  is  the  sum  of  the  original  operator 
A  along  the  diagonal  and  additional  terms  that  arise  due  to 
partial  derivatives  of  ^  on  the  right-hand  side.  Note  that  the 
representation  is  not  unique,  in  that 


^(0  = 


-ML 

dx‘^ 

COs{^x)  0^ 


0 

a" 

dx^ 


is  also  an  admissible  representation.  It  may  be  possible  to  use 
convex  combinations  of  admissible  representations  in  order  to 
elicit  better  global  convergence  properties  (since  each  repre¬ 
sentation  is  generally  nonlinear)  as  in  [5],  but  this  possibility 
is  not  explored  in  the  present  paper. 


We  prove  the  boundness  of  Cwp,  which  is  essential  to  ensure 
the  boundedness  of  the  control  signal  u{t)  in  (8).  This  result 
is  a  direct  generalization  of  Lemma  3. 


Lemma  7.  If  there  exist  constants  pi^  p2  >  0  such  that 
||w(t)||z  <  pi  and  ||wa^(t)||z  <  p2  for  some  t  >  0,  then 
there  exist  constants  kq  =  /^o(pi5  P2)  ond  ki  =  a:i(pi,  P2) 
such  that 


ma.x{\\wp{t)\\z,  \\wp,x{t)\\z)  <  Ko\\w{t)\\z  +  ki  (14) 

SO  that 

Ivpitp  =  \Cwp{t)f  <  2{Ko\\w{t)\\z  +  Ki)  (15) 

Proof.  To  prove  (11),  we  note  that  the  dynamics  of  Wp^x 
can  be  written  as 


Wp,xt  =  Awp^x  +  9{t,  X,  w) 

where 

df{t,x,w)  df{t,x,w) 

i,(..x,w)=  w. 

Applying  Lemma  1  and  taking  the  norm  of  both  sides  yields 
(14),  from  which  (15)  follows  as  a  direct  consequence  of 
applying  Agmon’s  inequality  (Lemma  2).  ■ 

We  state  the  main  result  and  allude  to  the  fact  that  its  proof 
is  identical  to  Theorem  1. 


Theorem  2.  Consider  the  system  (13).  Suppose  that  a  stable 
representation  ^(w)  exists  and  furthermore,  there  exists  an 
operator  Aq  which  gives  a  lower  bound  on  the  decay  rate  of 
trajectories  corresponding  to  ^(w).  Then  there  exists  a  semi¬ 
group  To{t)  corresponding  to  Aq  satisfying  the  growth  bound 
similar  to  Assumption  2.  Furthermore,  there  exists  a  BIBO  ( in 
the  sense  of  Coo)  stabilizing  controller  given  by  (5),  (6),  and 
(8)  provided  the  small  gain  condition  of  Theorem  I  is  satisfied 
{with  the  operator  T{f)  replaced  by  To(t)). 

Although  we  derived  the  above  results  for  a  class  of  systems 
of  the  form  w  =  Aw  +  /(t,  x,  w)  with  A  stable,  the  method 
can  be  readily  used  for  systems  of  the  form  w  =  AgW  F 
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/(t,x,w),  where  Ag  need  not  be  stable  and  could  even  be 
zero.  The  system  be  rewritten  as 

w  =  Aw  +  h{t,  X,  w),  h{t,  X,  w)  =  /(t,  X,  w)  +  AgW  —  Aw 

(16) 

and  Theorem  2  can  be  employed  directly.  Conditions  under 
which  a  particular  A  (as  compared  to  some  stable  A)  exists 
is  an  open  problem. 

V.  Insights  from  Finite  Dimensional  ODE  Systems 

In  this  section,  we  will  start  with  a  finite-dimensional 
linear  time-invariant  (LTI)  system  to  identify  conditions  under 
which  a  DPO  controller  would  lead  to  a  stable  closed-loop 
system.  Thereafter,  we  will  specialize  Theorem  1  to  finite 
dimensional  nonlinear  systems.  The  resulting  controller  and 
the  small  gain  stability  condition  will  be  similar  to  those 
obtained  in  adaptive  control  [7].  The  conditions  obtained  in 
this  section  would  be  important  while  using  DPO  in  a  practical 
setting,  where  hardware  implementation  and  verification  and 
validation  protocols  would  run  in  a  finite  dimensional  setting. 

A.  Application  to  Linear  Time -Invariant  Systems 
Consider  the  ODE  system 

w{t)  =  Aw{t)  -f  Bu{t),  y{t)  =  Cw{t)  +  Du{t)  (17) 

where  w  e  and  u^y  e  R.  The  system  is  assumed  to  be 
observable  and  controllable.  The  control  objective  is  to  ensure 
that  y{t)  tracks  a  reference  signal  r{t). 

The  particular  and  the  homogenous  halves  are  given  by 

Wp{t)  =  AmWp{t)  -  BKw{t),  yp{t)  =  Cwp{t) 

Wh(t)  =  AmWh(f)  +  Bu(t),  yh(t)  =  Cwh(t)  +  Du(t)  (18) 

where  Aj^  =  A  BK  is  Hurwitz.  Let  w  =  Wp  w^  —  w 
denote  the  error  between  the  observer  in  (18)  and  the  system 
in  (17).  The  dynamics  of  w  are  given  by  w{t)  =  Amw{t)  and 
y{t)  =  Cw{t).  Clearly,  w{t)  ^0  and  y{t)  -^0  exponentially. 

Let  G{s)  =  C{sl  -  Am)~^B  and  Gc{s)  =  G{s)  A  D. 
Then,  from  (18),  it  follows  that  yh{s)  =  Gc{s)u{s)  and 
yp{s)  =  —G{s)Kw{s).  A  control  signal  analogous  to  (8) 
is  given  by  u{s)  =  H{s){r{s)  -  yp{s)),  so  that  yh{s)  = 
Gc{s)H{s){r{s)-yp{s)).  We  want  yh{t)  to  track  r{t)-yp{t). 
The  ideal  choice  for  H{s)  would  then  be  1/Gc{s)  which  is, 
in  general,  neither  proper  nor  stable.  Therefore,  we  choose 
H{s)  so  that  i7(0)Gc(0)  =  1  to  guarantee  asymptotic  tracking 
and  impose  further  conditions  on  H{s)  to  ensure  a  stable 
closed-loop  system.  The  closed-loop  dynamics  are  obtained 
by  substituting  for  u{t)  in  (17),  whose  Laplace  transform  is 
given  by  sw{s)  =  Ax{s)-\-BH{s){r{s)-\-G{s)Kw{s)).  Since 
A  =  Am  —  BK,  we  get 

w{s)  =  {sI-AmRBK{l-G{s)H{s)))-^H{s)Br{s)  (19) 

Theorem  3.  The  closed-loop  system  (17)  and  (18)  is  expo¬ 
nentially  stable  if  and  only  there  exist  K  G  RA  and  a  stable 
transfer  function  H{s)  such  that:  (a)  A-\-BK  is  Hurwitz,  and 
(b)  all  poles  of  (4  +  {si  -  Am)~^BK{l  -  G{s)H{s)))-^ B 
lie  in  the  open  left  half  plane,  where  G{s)  =  G{sl  —  Am)~^B 
and  In  is  the  identity  matrix  of  size  n. 


Proof.  We  require  Am  to  be  Hurwitz  so  that  the  the  observation 
error  {w)  dynamics  are  stable.  Part  (b)  can  be  proved  by 
writing  the  closed-loop  system  in  the  form  (19).  ■ 

B.  Application  to  Finite  Dimensional  Nonlinear  Systems 
Next,  we  assume  that  the  system  is  of  the  form 

w  =  AmW  Fhu  A  f{w),  y  =  cw  (20) 

where  w  eR'^,  f  :R'^  ^  R'^,  and  u,  y  eR.  The  nonlinearity 
f{w)  is  generally  not  matched  to  u{t),  i.e.,  it  is  not  always 
possible  to  write  f{w)  =  bg{w)  for  some  nonlinear  g{w). 

Remark  1.  The  system  could  be  nonlinear  and  of  the  form 
w  =  fo{w)  +  bu,  in  which  case,  in  (17),  f{w)  =  fo{w)  — 
AmW.  Clearly,  Am  is  also  a  design  element,  constrained  by 
the  stability  and  performance  requirements. 

1 )  Perturbation  Observer  Design:  Let  Wp  and  Wh  denote 
the  states  of  the  particular  and  the  homogeneous  half,  and  let 
w  =  Wp  Wh.  Let  w  =  w  —  w.  The  dynamics  of  the  two 
halves  are  given  by 

Wp  =  AmWp  -h  f{w),  %  =  CWp  (21) 

Wh  =  AmWh  +  bu,  yh  =  cwh  (22) 

where  f{w)  is  the  predicted  value  of  f{w).  This  can  be  done 
using  a  wide  range  of  methods,  such  as  the  projection  operator 
and  we  retain  this  notation  for  simplicity. 

2)  Design  of  Control  Signal:  As  in  (8),  we  choose  a  control 
signal  u{t)  whose  Laplace  transform  is  given  by 

■u(s)  =  H{s){r{s)  -  yp{s))  (23) 

so  that  yhis)  =  c{sl  -  Am)~^bH{s){r{s)  -  yp{s}),  and  we 
prescribe  77(0)  =  —l/{cAmb). 

3 )  Closed-Loop  Stability:  In  order  to  show  stability  of  the 
closed-loop,  we  need  to  introduce  some  notation,  following 
that  introduced  earlier  for  the  general  system  represented  in 
the  operator  form. 

Assumption  5.  If  for  some  t  >  0,  \\w\\c^^t  <  Ph^  where  G 
IR+,  then  there  exist  constants  —  ^{ph)  —  ^i{Pb)y 

with  6o,  Si  G  M+  such  that  \\f{w)\\c^,t  <  So\\w\\c^,t  A  Si 
for  all  t  G  M. 

The  next  assumption  concerns  f{w)  in  (21).  Using  tech¬ 
niques  such  as  the  projection  operator  to  design  f{w),  we 
ensure  that  \\f{w)  —  f{w)\\oo  is  globally  uniformly  bounded. 

Assumption  6.  If  for  some  t  >  0,  <  Pb,  where 

pi)  G  then  there  exist  constants  uq  =  n{pi))  and  ni  = 
fs^i{pb)^  with  kq,  Ki  G  such  that  \\f{w)  —  f{w)\\c^,t  ^ 

Ko\\w\\c^)  +  Ki. 

Finally,  we  state  a  small-gain  design  condition  which  is 
critical  to  the  stability  proof  that  follows. 

Assumption  7.  There  exists  a  low  pass  filter  H{s)  and  a 
scalar  F  >  0  such  that 

1  -  ||(s/  -  A^)-^bH{s)\\cA\{sI  -  ^„^)-i£j|c||oo(/to  +  <^o) 
-  ||(s/-A„)-i|U,5o>r  (24) 
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(a)  Polynomial  design 


(b)  Finite  state  approximation  design 
Fig.  2.  DPO  illustrated  with  two  different  control  synthesis  methods. 

Define  a  constant  p  =  N/D,  where  the  numerator  N  is 
given  by 

N  =  ||(s/  -  A^)-HH{s)\\cMc^  (25) 

+  ^lUi  l|c||oo(<5l  +  f«l))  +  ||(sr  —  Am) 

and  the  denominator  D  by 

D  =  l-\\{sI-A^)-^UJo  (26) 

-  ||(s/  -  A^)-^bH{s)\\c,  ||(s/  -  Am)-'||A  ||c||oo((5o  +  Ko) 

We  state  the  main  results,  similar  to  Lemma  3  and  Theo¬ 
rem  1.  The  proofs  are  omitted  for  brevity. 

Lemma  8.  Consider  the  particular  half  of  the  observer  in  (21). 

II^IUoo  ^  P  —  ^ith  N  and  D  defined  in  (25)  and 

(26),  then  the  observer  state  Wp  satisfies  the  bound  < 

“  ^m)  ^||£i((^0  +  ^o)||'f^||£oo  +  (^1  +  ^l))- 

Theorem  4.  The  closed-loop  system  (20),  (21),  (22),  and  (23) 
is  BIBO  stable  and  the  state  w  of  the  closed-loop  system 
satisfies  ||w||£^  <  p. 

VI.  Simulations 
Consider  the  forced  wave  equation 

Quit,  x)  -  0.10txx{t^  x)  -  20yy{t,  x)  =  10006>(t,  x)  (27) 

pO.l 

Ox{t,  L  =  0.1)  =  0^  0{t,0)  =  u{t),  C0{t^x)=  /  0(t^x)dy 

Jo 

Note  that  the  dynamics  in  (28)  are  unstable.  The  DPO  is 
designed  as  follows: 

Op,tt  -  O.lOp^txx  -  ‘^Op^xx  =  F{t,x)  -p{0.10t{pfi)  H-  2§p{t)) 
Ktt  -  0.10h,txx  -  20h,xx  =  -p{0.10t{K  t)  +  20h{t))  (28) 

Op{t,  0)  =  0p^x{t,  0.1)  =  0,  0h{t,  0)  =  u{t),  0h,x{t,  0.1)  =  0, 

where  p  is  a  design  parameter.  The  controller  is  designed  using 
two  methods,  and  the  simulations  for  each  design  are  shown 
in  Fig.  2.  The  first  controller  is  designed  by  approximating 


the  homogeneous  half  of  the  observer  corresponding  to  (28) 
by  a  polynomial,  as  explained  in  [12].  The  second  controller, 
of  the  form  (23),  is  derived  for  a  finite  state  representation 
of  (28)  obtained  using  Galerkin’s  method.  In  both  cases,  as 
predicted,  the  closed-loop  system  is  stable  and  the  tracking 
error  is  almost  negligible.  The  closed-loop  system  is  stable 
when  p  G  [450,  550]  (approximately).  This  is  a  consequence 
of  having  chosen  the  right-hand  side  of  (28)  in  that  particular 
form.  The  polynomial  expansion-based  controller  presented 
in  this  paper  has  also  been  demonstrated  in  experiments  on 
the  bending  vibrations  of  a  beam  [11],  [12].  The  closed-loop 
stability  criterion  does  not  necessarily  inform  what  a  good 
choice  for  the  RHS  would  be,  and  an  optimum  design  of  the 
RHS  remains  an  open  problem. 

VII.  Conclusion 

In  this  paper,  we  presented  a  novel  control  method  for 
finite  as  well  infinite  dimensional  systems  based  on  a  two- 
stage  perturbation  observer.  The  control  signal  is  designed 
for  the  homogeneous  half  of  the  observer,  while  “perturbation 
prediction”  is  accomplished  by  putting  together  both  halves. 
The  robustness  of  the  controller  was  proved  using  the  small 
gain  theorem.  In  particular,  we  showed  that  the  £oo-nofin  of 
the  system  state  (over  time  for  finite  dimensional  systems,  and 
over  space  as  well  as  time  for  PDFs)  as  well  as  the  control 
inputs  is  uniformly  bounded.  Simulations  were  performed  to 
demonstrate  the  effectiveness  of  the  controller. 
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1.  Introduction 


For  the  last  decades,  scientists  are  highly  interested  in  designing  biomimicry  devices  performing 
as  artificial  muscles.  Compliant  and  light\A/eight  robotic  arms  and  grippers  for  tight  spaces  as  well 
as  energy-efficient  actuators  for  biomimetic  locomotion  have  drawn  much  attention  for  study. 
To  develop  a  practical  biomimic  actuator  with  good  performance,  choosing  suitable  materials  is 
one  of  the  most  important  and  fundamental  issues.  Smart  materials  that  can  be  triggered  to 
deform  through  certain  stimulations  are  what  we  desire.  Besides,  muscle-like  features,  including 
light,  soft,  and  short  responding  time,  are  also  expected. 

In  literatures,  numerous  of  smart  materials  being  proposed.  From  these  materials,  electroactive 
polymers  (EAPs)  with  lightweight  and  nice  compliancy  are  very  suitable  to  serve  as  candidates 
for  biomimic  actuator  applications.  EAPs  are  polymers  that  are  cable  of  inducing  deformation 
under  electrical  stimulation.  The  human-muscle  like  features,  such  as  lightweight,  considerable 
large  displacement,  and  acceptable  response  time  provide  EAPs  high  potential  of  artificial  muscle 
applications. 

A  group  of  the  EAP  family,  dielectric  elastomers  (DEs),  has  especially  drawn  our  attention  due  to 
their  high  strain,  comparably  short  response  time,  low  cost,  and  high  electromechanical  coupling 
efficiency  [1-5].  They  are  often  used  to  fabricate  actuators,  so  called  dielectric  elastomer 
actuators  (DEAs).  Basically,  DEAs  are  made  of  incompressible  soft  dielectric  elastomer 
membranes  sandwiched  between  compliant  electrode  layers  to  form  dynamic  capacitors.  When 
electric  field  is  applied  across  the  electrodes,  the  columbic  force  generates  Maxwell  stress  [6] 
that  attracts  other  electrodes  together  and  squeezes  the  sandwiched  incompressible  dielectric 
elastomer  layer.  As  a  result,  the  in-plane  expansion  of  DEA  can  be  observed. 

Many  different  designs  of  soft  actuator  have  been  developed  with  the  usage  of  dielectric 
elastomers.  Two  major  ideas  were  presented  in  the  past:  a)  actuators  with  hard  frame  for 
supporting  prestretch,  and  b)  using  prestretch  as  stored  energy  for  the  recovering  of  bi-stable 
structures.  For  an  example  of  hard  frame  supporting  DEA,  a  robotic-arm  design  that  S.  Dubowsky 
et  al.  published  in  2006  [7]  uses  frames  with  deformable  joints  to  support  prestretch  and,  in  the 
meantime,  allowing  dielectric  elastomers  to  expand  while  they  electric  field  is  applied.  This 
structure  makes  a  semi-freestanding  actuation  that  frame-included  movement  is  presented 
instead  of  in-frame-only  displacement.  Several  cells  of  actuators  can  be  combined  together  for 
multi-degree-of-freedom  motion.  Similarly,  by  combing  DEAs  with  mechanical  structures,  P. 
Lochmatter  et  al.  [8]  presented  a  shell  bending  actuator  with  a  stack  of  segments.  The  shell-like 
structure  has  two  parallel  planar  actuators  and  joint  mechanical  parts  in  between.  Bending  can 
be  observed  once  one  side  of  DEAs  is  stimulated  and  expanded.  Also,  series  of  cells  can  be 
connected  together  to  induce  larger  bending  motion. 


On  the  other  hand,  one  of  the  bi-stable  structure  actuator  designs  named  dielectric  elastomer 
minimum  energy  system  (DEMES)  were  presented  by  G.  Kofod  [9].  With  prestretched  elastomers 
attached  on  deformable  frames,  actuators  may  stay  in  relaxation  phase  where  elastomers  are 
crumbled  with  no  electric  field  applied.  And  when  the  actuators  are  activated,  they  may  deform 
into  the  other  stable  phase  where  the  energy  from  applied  electric  field  overcomes  the  stored 
energy  provided  by  prestretch  and  makes  elastomer  flat.  Later,  Petralia  and  Wood  [10]  further 
presented  their  work  of  a  curving  chain  actuator  by  combining  several  paper-frame-based  DEMES 
cells. 

In  literatures  [11,12],  prestretch  is  a  method  showing  great  improvement  of  in-plane  actuation 
strain  of  dielectric  elastomer  actuators.  Dielectric  elastomer  actuators  with  prestretch  can 
achieve  more  than  300%  of  strain.  In  all  of  these  designs,  prestretch  is  taking  parts  in  the 
fabrication  that  benefits  actuators  having  large  deformation.  However,  the  additional  supporting 
frames  cause  operational  setbacks.  For  example,  weight  and  space  that  the  prestretch-support 
structures  provide  are  much  more  than  dielectric  elastomer  itself  [13,14].  In  addition,  non- 
uniform  prestretch  and  stress  relaxation  may  affect  subsequent  actuation  [15]  and  cause  local 
strain  division.  From  our  previous  experimental  observation  of  planar  DEAs  [16],  actuation  strain 
will  be  influenced  by  non-uniform  prestretch  and  may  be  controlled  by  the  biaxial  prestretch 
ratio.  Larger  prestretch  tends  to  induce  lower  actuation  strain,  and,  vice  versa,  lower  prestretch 
tends  induce  higher  actuation  strain. 

To  minimize  the  inconveniences  come  with  prestretch,  we  present  another  design  of  dielectric 
elastomer  actuator  fabrication  with  no  requirement  of  prestretch  and  no  additional  equipment 
for  supporting  prestretch.  In  this  study,  freestanding  design  of  DEAs  that  increases  the  energy 
mass  ratio  and  the  flexibility  of  actuators  will  be  revealed.  Our  idea  is  using  partial  surface 
reinforcement  on  thin  planar  dielectric  elastomer  actuators  to  induce  out-of-plane  deformation. 
When  electric  field  is  applied,  planar  actuators  tend  to  perform  in-plane  expansion  due  to  the 
Maxwell  stress  while  stiffeners  added  on  the  surface  act  as  surface  reinforcements  and  constrain 
expansion  of  designed  direction.  They  induce  out-of-plane  deformation  and  guide  the  bending 
motions.  By  applying  different  configurations  of  stiffeners  on  different  positions,  multiple  choices 
of  three  dimensional  deformations  with  different  directions  and  degrees  of  bending  and  twisting 
are  expectable.  Since  the  planar  actuator  is  based  on  thin  dielectric  elastomer,  a  considerable 
out-of-plane  deformation  doesn't  require  large  in-plane  expansion  that  relies  on  the  help  of 
prestretch. 

In  Suo's  work  [17],  they  showed  experimental  results  and  analytical  analysis  of  applying  fiber 
reinforcement  on  the  surface  of  dielectric  elastomer  actuator.  When  electric  field  is  applied,  the 
expansion  of  elastomer  is  constrained  on  the  direction  of  fibers  and  the  actuators  perform 
unidirectional  actuation.  Using  similar  idea,  stiffeners  in  our  work  with  different  are  attached  on 


one  surface  of  a  planar  actuator  partially  and  constraining  the  expansion  on  the  surface.  This  will 
cause  the  entire  planar  actuator  deforms  out-of-plane  significantly. 

In  this  paper,  we  will  show  the  fabrication  of  planar  actuators.  Analytical  analysis  and  FEM 
simulation  will  also  be  provided  to  investigate  the  role  of  stiffener  which  can  be  used  to  change 
the  force-stroke  characteristics  of  actuators. 


2.  Fabrication  Procedure 

The  fabricated  device  has  a  square  shape  of  25x25mm.  It  contains  elastomer  layers  (3M  VHB 
F9460PC  tapes)  and  electrode  layers  (carbon  black  powder,  Super  C65,  TIMCAL  Inc.,  USA).  The 
elastomer  layers  include  active  area  in  the  middle  with  22x22mm  and  inactive  border  of  1.5mm 
around  the  edge  for  sealing  and  preventing  circuit  shorting  (Figure  la).  A  window  mask  is  used 
to  define  the  electrode  area.  Carbon  black  powder  is  uniformly  brushed  over  the  window  mask. 
Narrow  strips  of  Aluminum  foils  are  attached  to  the  edge  of  electrode  layer  to  form  external 
terminals  of  the  actuator.  Finally,  another  VFIB  tape  is  applied  on  the  stack  as  a  cover.  The 
sequence  is  repeated  for  multiple  stack  modules.  The  aluminum  terminal  locations  are  alternated 
to  separate  positive  and  negative  electrode  when  fabricating  multi-layer  structure.  As  an  example, 
sketches  of  1  (which  is  the  minimum  stack  of  basic  DEA  design),  2,  and  3-cell  structure  are  shown 
in  Figure  2.  After  this,  the  stiffener,  3M  Magic  Scotch  tapes,  are  cut  into  long-narrow  shape  with 
3mm  width  and  attached  on  the  surface  of  the  DEA  laminates  stack  (Figure  lb).  Different  shapes 
and  configurations  of  stiffeners  are  able  to  apply  for  more  complicated  deformation. 


aluminum  foil 

sealing  area 
(width  1,5mm) 

Figure  1.  Sketches  of  DEA  samples,  (a)  Top  view  of  unit-cell  DEA  structure  and  its  measurement 
with  (b)  single  strip  configuration  stiffener  in  the  middle  on  the  top  surface. 


Figure  2.  Sketches  of  DBA  samples  fabrication  with  side  view  of  different  stack  configurations,  (a) 
1-cell  (b)  2-cell  (c)  3-cell. 


3.  Experiment 

In  the  experiments,  DEA  samples  were  hung  from  the  electrodes  terminal  and  connected  to  the 
high  voltage  circuit  terminal  and  applied  maximum  voltages  range  of  0  -  2.24kV  by  utilizing  a  DC- 
DC  voltage  converter  (Q-80,  EMCO  Inc.).  The  converter  has  a  high  DC  voltage  linear  amplifier 
range  of  0  -  8kV  output  for  0  -  5V  input.  The  actuator  deformations  were  captured  by  high- 
resolution  CCD  camera  (2448x2048  pixel,  Grasshopper,  Point  Grey  Inc.)  with  an  in-situ  image 
system  setup  for  out-of-plane  displacement  and  curvature  analysis. 

In  this  section,  experimental  results  of  1,  2,  and  3-cell  actuators  with  1-strip  configuration  of 
stiffener  are  shown  below. 
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Figure  3  shows  the  bending  deformations  of  the  actuators  with  the  increment  of  amount  of 
stacking  cell  from  the  side  view  under  applied  nominal  electric  field  of  0  -  22.4MV/m.  When  the 
input  nominal  electric  field  is  increasing  as  well  as  when  more  cells  are  laminated,  larger  bending 
effects  were  observed  consistently  in  three  different  actuators.  Planar  actuators  are  fabricated 
with  DBA  laminates  which  contains  multi  stack  of  dielectric  elastomer  actuators.  By  increasing 
the  number  of  laminated  cells,  higher  energy  is  provided  that  may  be  able  to  generate  larger 
deformation. 

Curvatures  were  calculated  from  the  captured  images  in  order  to  define  the  deformation  of 
actuators.  Partial  circles  were  assumed  for  every  bending  motion  to  ease  the  measurement  and 
calculation.  The  results  are  shown  in  Figure  6. 
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Figure  3.  Experimental  results  of  deformation  sequence  of  1,  2,  3-cell  DBA  with  1-strip 
configuration  of  stiffener  attached  under  applied  electric  field  of  0-22.4MV/m. 


In  addition,  surface  deformation  analysis  is  also  performed  by  using  digital  image  correlation 
technique.  The  surfaces  of  planar  actuators  were  patterned  with  alumina  powders.  Reference 
and  deformed  images  were  captured  before  and  after  actuation.  Images  were  analyzed  by  using 
commercial  software,  Vic2D  (Correlated  Solutions,  Inc.)  to  obtain  and  evaluate  the  in-plane  finite 
strain  components.  The  results  of  surface  deformation  analysis  will  be  used  to  compare  with 
finite  element  analysis  in  order  to  investigate  the  role  of  stiffeners. 


4.  Analytical  analysis 

To  understand  the  behavior  of  this  actuator  design,  we  utilized  the  Timoshenko's  analysis  of  bi¬ 
metal  thermostats  [18]  and  introducing  the  Maxwell  stress  instead  of  the  thermal  effect  to  arrive 
at  the  general  deformation  representation  under  applied  electric  field.  The  derivation  of  the 
analysis  is  illustrated  as  follows. 

A  free-body  diagram  of  stiffener  and  EAR  bi-layer  structure  can  be  simply  drawn  in  Figure  4.  When 
the  actuator  is  activated,  total  strain  of  each  layer  is  a  summation  of  strain  that  caused  by 
electrical,  axial,  and  bending  load.  It  can  be  written  in  equation  (1)  described  with  material 
modulus,  E  and  geometry  parameters  of  layer  /. 
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Figure  4.  Free-body  diagram  of  stiffener  and  EAR  bi-layer  structure  represents  planar  actuator 
with  surface  reinforcement.  In  this  figure,  a  is  the  thickness,  w  is  the  width,  M  is  the  bending 
moment,  and  R  is  the  axial  load. 


First  of  all,  to  define  electrical  expansion  of  EAR,  we  apply  Maxwell  stress  equations  shown  in 
equation  (2)  and  Flooke's  Law  assuming  in-plane  deformation  is  small  and  within  the  elastic 


deformation  zone.  Electrical  strain  is  obtained  and  written  as  a  function  of  applied  nominal 
electric  field  in  equation  (3). 
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From  the  definition  bending  moment  in  beam  theory,  m  =  —  ,  where  p  \s  the  radius  of  the 

P 

bended  beam.  Apply  this  into  equation  (1)  and  the  bending  term  of  strain  can  be  derived  and 
written  as  in  equation  (4). 
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Since  the  axial  load  =  W2P2  =  P  ■  Besides,  the  total  bending  moment  can  be  represented  by 
axial  load  multiple  by  hale  of  the  total  thickness  as  well  as  the  summation  of  bending  moments 

on  each  layer,  =  +M, .  Therefore,  axial  term  of  strain  can  be  written  as  in  equation 
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In  addition,  because  of  the  continuity,  total  strain  at  the  interface  of  layer  1  and  2  should  be  the 
same.  Therefore,  the  strain  condition  after  actuation  can  be  described  by  equation  (6). 
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Let  A  =  +  ^2  /  actuation  curvature  k  ,  which  is  equal  to  inverse  of  radius  p,  of  this  system  can 

be  written  as  a  function  of  material  modulus  and  geometry  parameters  in  equation  (7). 
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5.  Computational  simulation 

Finite  element  method  is  utilized  for  the  computational  simulations.  User  Material  subroutine  is 
applied  with  commercial  finite  element  software  ABAQUS.  Incompressible  Neo-Hookean  model 
is  used  to  represent  the  elastomer.  From  literature  [19],  we  know  that,  elastic  energy  and 
electrostatic  energy  of  dielectric  elastomer  actuator  can  be  written  as  a  function  of  stretch  1 
and  applied  nominal  electric  field  Oq  as  equation  (8)  and  (9). 
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Therefore,  the  total  free  energy  can  be  written  in  equation  (10) 
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Here,  nominal  electric  field  is  defined  as  input  voltages  used  in  experiments  divided  by  original 
thickness  between  electrode  layers  of  unit  cell  DBA  (0.1mm). 

As  a  result,  the  two  required  material  properties  in  our  finite  element  simulation  are  shear 
modulus  and  dielectric  constant.  The  dielectric  constant  of  VHB  9460PC  tapes  is  used  for  3.21  [20] 
and  the  shear  modulus  is  found  to  be  0.03  MPa  from  a  uniaxial  tensile  tests  we  performed  shown 
in  Figure  5.  Assuming  our  interest  deformation  is  small  and  in  the  elastic  behavior  region.  Young's 
modulus  is  calculated  to  be  0.1  MPa  from  the  strain-stress  curve.  Shear  modulus  is  calculated 
and  obtained  to  be  0.03  MPa  by  applying  Hooke's  Law  where  0.5  for  Poisson's  ratio  is  employed 
assuming  incompressible  material. 
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Figure  5.  Stain-stress  curve  of  VHB  F9460PC  tapes  from  uniaxial  tensile  test. 

In  the  computational  simulation,  3D  geometric  half  model  for  cantilever  laminate  was  used  for 
overall  structure.  Thickness  of  electrode  layers  is  neglected.  For  multiple-cell  actuator  cases,  all 
active  layers  were  combined  into  one  single  layer.  A  Hybrid  20-node  quadratic  brick  element 
employed  for  each  individual  layer  linked  with  surface-based  tie  constraint.  One  edge  of  the 
sample  was  anchored  and  other  edges  were  left  free-to-move.  Curvatures  are  used  to  represent 
actuator  deformation  in  simulation.  They  were  assumed  as  partial  circles  and  calculated  by  using 
the  center  in-plane  and  out-of-plane  displacement.  Figure  6  shows  the  comparison  of 
computational  and  experimental  results. 

From  finite  element  analysis,  the  profile  of  in-plane  strain  field  is  obtained  to  compare  with 
experimental  observation  from  surface  displacement  analysis.  Results  from  both  analyses  (Figure 
7)  show  a  band  of  region  that  much  higher  strain  occurs.  The  high  strain  region  represent  lateral 
strain  perpendicular  to  the  stiffener  directon  concerntrates  near  the  boundary  of  stiffeners.  The 
concerntrated  strain  is  caused  by  stiffener  constraints  to  the  elastomer  expansion.  This  shows 
the  nonliearity  of  in-plan  deformation  of  this  type  of  actuator  which  will  affect  on  our  analitical 
analysis.  More  detail  will  be  discussed  in  the  next  section. 
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Figure  6.  Curvature  comparison  of  experimental  and  computational  results  of  1,  2,  3-cell 
actuators  with  1-strip  configuration.  Solid  lines  represent  FEM  results  and  square  marks 
represent  experimental  results. 
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Figure  7.  Surface  displacement  analysis  (a)  and  finite  element  (b)  results  of  concentrated  strain 
around  stiffener  under  applied  nominal  electric  field  of  20800  V/m. 


6.  Discussion 


For  this  special  design,  two  cover  layers  are  attached  and  their  effects  to  the  electrical  strain  have 
to  be  considered.  Therefore,  a  parameter  p  is  requred  in  the  developed  analytical  solution  to 
represent  how  the  additonal  symmetric  cover  layers  influence  the  in-plane  expansion,  p  is 
defined  as  the  thichness  of  active  layers  devided  by  the  total  thickness  of  actuators.  For  example, 
for  case  of  1,  2,  and  3-cell  actuators,  p  is  1/2,  2/3,  and  3/4  respectively. 

Comparing  our  analytical  solution  with  experimental  data,  a  fitting  parameter  in  the  analytical 
solution  is  found  in  order  to  represent  all  of  the  nonliearities  including  inhomogeous  deformation 
and  nonliear  electromechanical  coupling  effect.  In  this  case,  the  fitting  parameter  is  determined 
to  be  0.55.  Finally,  the  analytical  solution  can  be  written  in  equation  (11).  Actuation  curvaures 
are  reprensented  as  a  function  of  input  nominal  electric  field. 


0.55/? 

^  1 

^ - 

-{E,I,+EJ,) 


1 


-I-- 


(11) 


Using  this  developed  analytical  solution,  the  deformation  sequences  under  serial  applied  electric 
field  from  experiments  are  compared  and  plotted  in  Figure  8. 
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Figure  8.  Curvatures  from  deformation  sequences  versus  applied  electric  field  \A/here  the  dotted 
lines  are  analytical  solutions  and  square  marks  are  experimental  results.  The  color  blue,  red,  and 
black  represent  1,  2,  and  3  cell  actuators  respectively. 


Force-stroke  characteristic  of  can  be  obtained  from  finite  element  analysis.  In  the  computational 
simulation,  we  first  constrain  the  out-of-plane  displacement  of  the  far-end  edge  opposite  to  the 
anchored  edge  and  then  applied  electric  field  on  the  actuators  to  certain  level.  After  the  electric 
field  reaches  certain  level,  the  edge  constraint  is  removed  and  replace  with  uniform  out-of-plane 
displacement  applying  on  tide-up  edge.  Reaction  forces  on  every  node  along  the  edge  are 
recorded  and  summed  up  in  every  interested  displacement  step.  The  results  of  the  force-stroke 
characteristic  of  a  3-cell  actuator  stiffener  surface  reinforcement  under  different  applied  nominal 
electric  fields  are  shown  in  Figure  9. 

In  Figure  9,  the  solid  lines  show  the  force-stroke  characteristics  of  a  3-cell-l-stiffener  actuator 
are  controlled  by  applied  electric  field.  With  increasing  applied  electric  field  to  higher  level, 
larger  reaction  force  is  observed  at  zero  displacement  and  larger  displacement  is  observed  at 
zero  reaction  force.  However,  for  a  free-to-move  actuator,  nodes  at  edge  are  not  tied  up.  This 
leads  to  deviation  of  reaction  force  calculated  from  the  summation  of  each  node  and  the 
experimental  data  converted  from  actuation  curvatures  which  are  represented  in  solid  circles  in 
the  plot. 

Figure  9.  Force-stroke  characteristics  of  a  3-cell-l-stiffener  actuator  under  applied  electric  field 
of  17600  (green),  19200  (red),  and  20800  V/mm  (blue). 

7.  Conclusion 

In  this  paper,  we  shared  an  idea  of  freestanding  planar  actuators  based  on  dielectric  elastomers 
capable  of  out-of-plane  deformation  by  using  partial  stiffener  reinforcement  on  the  surfaces. 
With  our  concept,  no  additional  devices  to  support  prestretch  are  required.  Deformation  of  the 
actuator  can  be  designed  is  by  using  different  configurations  of  surface  stiffener  reinforcements. 
Experiments  provided  the  actuators  performances.  In  addition,  surface  displacement  analyses 
were  employed  in  order  to  assist  understanding  how  stiffeners  induce  deformation  to  planar 
actuators.  An  effective  band  around  stiffener  constraints  was  discovered  which  leads  to 
inhomogeneous  deformation  in  this  case.  The  analytical  analysis  was  provided  for  future  control 
strategies.  Finite  element  framework  was  also  developed  which  could  be  used  to  investigate  the 
role  of  stiffeners  and  force-stroke  characteristics  of  actuators.  It  could  also  be  employed  to  study 
more  complex  3D  motion  based  on  this  actuator  design  in  our  future  works. 
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Introduction 


Dielectric  elastomers  (DEs)  are  a  type  of  electroactive  polymers  that  has  been  studied  for  many 
years  in  the  booming  field  of  biomimicry  devices  development.  The  high  strain,  comparably  short 
response  time,  low  cost,  and  high  electromechanical  coupling  efficiency  []  make  them  a  strong 
candidate  of  material  selection  as  artificial  muscles.  They  are  often  used  to  fabricate  actuators, 
so  called  dielectric  elastomer  actuators  (DBAs).  DBAs  are  made  of  soft  dielectric  elastomer 
membranes  sandwiched  by  compliant  electrode  layers  to  form  dynamic  capacitors.  The  columbic 
force  generates  Maxwell  stress  []  when  electric  field  is  applied  across  the  thickness  direction. 
Maxwell  stress  attracts  both  electrode  layers  together  and  compress  the  middle  incompressible 
elastomer  which  will  introduce  the  in-plane  expansion  of  DBAs. 

The  idea  of  fabricate  planar  actuators  based  on  dielectric  elastomer  has  been  discussed  in 
electroactive  polymer  society  for  years  [].  Due  to  the  mechanism  of  dielectric  elastomer  actuators, 
motion  of  in-plane  expansion  is  very  direct  and  easy  to  achieve.  However,  in  order  to  broaden 
the  application  possibilities  and  categories,  it  is  important  to  study  structures  that  can  induce  3- 
dimentional  motion.  Several  different  designs  were  presented  in  the  literatures. 

Since  it  has  been  shown  that  prestretch  is  a  method  capable  of  greatly  improving  in-plane 
actuation  strain  of  dielectric  elastomer  actuators  [],  a  branch  of  planar  dielectric  elastomer 
actuators  design  is  based  on  prestretched  membrane.  P  Lochmatter  and  G  Kovacs  demonstrated 
an  idea  of  shell-like  actuators  []  []  using  two  parallel  prestretched  dielectric  elastomer  actuators 
connected  by  hinge  structures.  When  the  actuator  on  one  side  of  the  shell  is  activated,  shell  is 
capable  of  bending.  The  core  structure  of  hinge  can  also  be  constructed  in  different  ways,  such 
as  an  array  of  spheres.  This  may  help  the  shell-like  actuator  present  biaxial  deformations  for  more 
complicated  motion.  In  their  later  studies  [],  the  actuator  can  be  further  improved  by  several 
adjustments,  for  example,  multi  layers  dielectric  elastomers  and  reinforced  carbon  rods  which 
support  the  prestretch.  Thanks  to  the  prestretch,  these  bending  devices  showed  results  of  large 
out-of-plane  deformation  and,  more  impressively,  outstanding  blocking  force  from  the  bending 
motions.  Besides,  because  of  the  prestretch,  higher  input  electric  field  can  be  applied  and,  as  a 
result,  actuators  demonstrate  even  better  bending  characteristics.  However,  the  additional 
supports  also  cause  unwanted  setbacks  that  conflicts  the  goal  of  developing  light  weight  and 
flexible  actuators  when  employing  electroactive  polymers  at  the  first  place.  For  example,  support 
structures  for  prestretch  occupy  much  more  weight  and  space  than  elastomer  itself  []  and  limited 
the  degree-of-freedom  of  deformations  by  the  supporting  and  connecting  structures.  This  branch 
of  design  is  featuring  high  value  performances  of  deformation  with  a  trade-off  of  losing  flexibility 
and  energy-to-mass/space  ratio  by  using  additional  rigid  prestretch-support  frames. 

Another  approach  to  induce  out-of-plane  deformation  from  planar  actuators  is  named  dielectric 
elastomer  minimum  energy  system  (DEMES).  The  idea  was  first  presented  by  G  Kofod  et  al.  [] 


which  demonstrated  a  clever  strategy  reversibly  utilizing  prestretch.  Dielectric  elastomers  are 
prestretched  and  attached  on  supporting  structures  which  are  inextensible  but  pliant.  Usually 
thin  membranes  of  plastic  slides  or  papers  are  employed.  Without  additional  energy  input,  the 
elastic  energy  used  to  prestretch  DE  films  will  transfer  to  the  outer  frame  and  causes  the  entire 
system  to  buckle.  Once  electric  field  is  applied  to  the  system,  the  expansion  of  elastomer  will 
flatten  the  system  and  reverse  the  buckling.  The  system  motions  are  defined  by  the  frame 
structure.  Different  configurations  of  frames  may  induce  different  shapes  of  buckling.  Due  to  the 
usage  of  prestretch,  deformation  and  output  blocking  force  of  the  system  is  comparably  large  []. 
G  Kofod  et  al.  []  showed  an  application  of  griping  actuator  in  their  study  demonstrated  the 
significant  force  that  this  system  may  provide.  Besides,  the  pliant  frames  give  the  system  higher 
flexibility  comparing  to  rigid  frames  and  extend  broader  application.  Despite  the  beauty  features 
of  DEMES,  a  disadvantage  appears  in  deformation  shape  prediction  with  difficulty  and  inaccuracy. 
M  T  Petralia  and  R  J  Wood  presented  a  structure  for  highly-deformable  soft  robotic  system  based 
on  this  strategy.  In  the  pictures  from  their  study,  we  can  observe  noticeable  differences  of 
deformation  even  with  the  same  support-frame  shape. 

Other  than  the  above  designs,  another  approach  is  the  unimorph  DEAs.  In  the  early  study  of 
dielectric  elastomer  actuators  [],  a  simple  idea  to  provide  unimorph  bending  effect  by  attaching 
extendable  DE  with  stiff  substrate  was  proposed.  In  2011,  O  A  Araromi  et  al.  []  []  introduced  a 
comprehensive  study  cover  from  an  unique  spay-deposited  fabrication  to  actuators 
characterization.  In  this  their  design,  multiple  layers  of  DE  film  were  spay  deposited  on  a 
substrate  membrane  with  compliant  electrode  in  between  each  layer  and  no  employment  of 
prestretch.  The  stiffer  substrate  constrained  one  side  of  the  planar  actuator  expansion  and 
induced  overall  bending.  This  bending  device  showed  considerable  deflection  but  less  blocking 
force.  Nevertheless,  without  the  requirement  of  prestretch  supporting  frames,  the  actuator  is 
lighter  and  more  flexible  and  still  performing  noticeable  deformation. 

To  not  only  keep  the  existed  advantages  of  unimorph  DEAs  such  as  light  weight  and  large 
deflection  but  also  further  improve  the  flexibility  of  actuator  motions,  we  introduced  a  strategy 
to  fabricate  planar  dielectric  elastomer  actuators  capable  of  out-of-plane  deformation  with  no 
requirement  of  prestretch  [].  Instead  of  applying  fully  covered  stiffer  substrate  to  urge  bending, 
partially  reinforcement  of  stiffener  segments  were  used  in  our  design.  Small  segments  of  stiffener 
reinforcement  are  able  to  provide  even  more  complicated  3-dimentional  deformation  to  a  planar 
actuator  and  at,  the  meantime,  featuring  light  weight  and  high  energy-to-mass  ratio.  In  this  study, 
we  will  provide  deeper  understanding  of  the  role  of  stiffeners  and  their  interactions  with 
dielectric  elastomers.  By  varying  stiffener  configurations,  we  are  able  to  change  the 
characteristics  of  actuators,  such  as  deflection  and  reaction  force,  while  keeping  bending  stiffness 
the  same.  This  enlarges  the  complexity  of  surface  motion  to  design  planar  actuators. 


Frabrication  and  Experiment 


To  study  the  role  of  stiffener,  we  planed  to  fix  the  geometry  parameter  of  bending  stiffeness  and 
vary  perodicity  of  stiffeners.  Two  cases  were  compared:  (a)  sigle  stiffener  and  (b)  split  siffener 
while  maintainting  the  totale  width.  Following  the  detail  of  fabrication  procedure  in  our  privious 
work  [],  two  3-cell  planar  actuators  were  fabricated.  A  3-cell  planar  actuator  is  composed  by 
stacking  three  unit  cells  of  actuator.  Each  unit  cell  includes  two  layers  of  VFIB  9460PC  tape  (50pm 
thick  each)  in  the  middle  and  sandwitched  by  compliant  electrode  layers  where  carbon  black 
particles  (Super  C65,  TIMCAL  Inc.,  USA)  were  employed.  Two  additional  VHB  tapes  were  used  to 
cover  up  the  stacked  actuator  for  protection  that  makes  a  total  thickness  of  0.4mm  for  a  3-cell 
planar  actuator  (thickness  of  electode  layers  were  neglected).  One  of  the  planar  actuators  was 
attached  with  one  3mm-wide  stiffener  (3M  Magic  scotch  tapes)  on  the  surface  and  the  other 
actuator  was  attached  with  three  Imm-wide  stiffeners  equally  seperated  as  shown  in  Figure  lb 
and  Ic.  Both  actuators  were  connected  to  the  high  voltage  circuit  with  applied  voltages  range  of 
0  -  2.24kV  by  employing  a  DC-DC  voltage  converter  (Q-80,  EMCO  Inc.).  High-resolution  CCD 
camera  (2448x2048  pixel.  Grasshopper,  Point  Grey  Inc.)  with  an  in-situ  image  system  setup  was 
utilized  to  capture  the  side  view  of  actuator  deformations.  The  deformation  sequence  of  3-cell 
DEA  with  stiffener  configuration  of  1  and  3  segments  under  applied  nominal  electric  field,  ,  of 

0-22.4kV/mm  is  shown  in  Figure  2.  Besides,  actuation  curvatures  were  measured  from  the 
captured  images  and  shown  in  square  marks  in  Figure  3. 
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Figure  1.  Schematic  diagram  of  (a)  side  view  of  3-cell  planar  actuator  and  two  DEA  samples  with 
(b)  singe  stiffener  and  (c)  triple  split  stiffeners  attached. 
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Figure  2.  Experimental  results  of  deformation  sequence  of  3-cell  DEA  with  stiffener  configuration 
of  1  and  3  segments  attached  under  applied  electric  field  of  0-22.4MV/m. 
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Figure  3.  Actuation  curvatures  under  applied  electric  field  of  0-22.4kV/mm  of  1  (black)  and  3  (red) 
segments  of  stiffeners  attached  on  a  3-cell  DEA.  Square  marks  and  dotted  lines  represent  results 
from  experiments  and  finite  element  analysis  respectively. 


Computational  simulation 


Utilizing  the  framework  we  developed  that  simulates  the  coupling  electromechanical  behavior  of 
dielectric  elastomer  actuators  [],  we  performed  in  computational  works  of  both  single  and 
splitting  stiffeners  cases.  In  the  simulation,  overall  structure  was  based  on  half  model  of  3D 
geometric  cantilever  laminate.  All  active  layers  were  combined  into  one  single  layer  and,  in  the 
meantime,  thickness  of  electrode  layers  was  neglected.  A  Hybrid  20-node  quadratic  brick 
element  was  employed  for  each  individual  layer  linked  with  surface-based  tie  constraint.  One 
edge  of  the  sample  was  anchored  and  other  edges  were  left  free-to-move.  Curvatures  were  used 
to  represent  actuation  motion.  Curvatures  were  assumed  as  partial  circles  and  calculated  based 
on  the  in-plane  and  out-of-plane  displacements  of  the  central  point  at  the  end  free  edge.  The 
deformation  results  of  both  single  and  three  stiffeners  cases  of  actuator  are  shown  in  Figure  3 
with  respect  of  serial  input  nominal  electric  field  up  to  22.4  kV/mm.  The  side  view  of  deformation 
sequence  in  simulation  is  shown  in  Figure  4.  Strain  field  of  the  actuators  after  electric  force 
stimulated  can  also  be  obtained  and  shown  in  Figure  5. 


Figure  4.  Deformation  sequence  in  finite  element  simulation  of  3-cell  planar  actuator  with  1  and 
3  segments  of  stiffeners  under  nominal  electric  field  of  0-22.4  kV/mm. 
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Figure  5.  Strain  field  of  3-cell  planar  actuator  with  1  segment  (a)  and  3  segments  (b)  of  stiffener 
under  nominal  electric  field  of  22.4  kV/mm. 

Discussion 

From  the  experiments,  we  observed  that  the  actuation  curvatures  are  controled  by  the 
configuration  of  attached  stiffeners.  In  the  two  cases  we  experimentally  performed,  the  total 
bending  stiffness  was  controled  to  be  the  same  by  maintaining  the  total  width  of  stiffeners. 
Flowever,  different  characteristics  were  revealed  in  two  actuators.  With  increasing  precocity  of 
stiffeners,  larger  curvatures  can  be  achieved.  To  investigate  this  phenomenon,  we  go  through 
finite  element  simulation  from  several  aspects.  First  of  all,  if  we  look  at  the  general  deformation 
shapes  shown  in  Figure  6,  bending  is  not  uniformly.  Buckling  occurred  between  the  stiffener 


edges  and  the  end-edge  of  actuator.  For  a  bending  device,  this  kind  of  buckling  is  not  what  we 
appreciate.  It  wastes  energy  on  unwanted  direction  (y-axis)  deformation.  As  a  matter  of  fact, 
increasing  precocity  of  stiffener  is  able  to  limit  the  unwanted  deformation.  And  doing  by  this, 
actuation  closer  to  homogeneous  bending  will  be  observed.  And  because  of  deformation  occurs 


on  the  bending  direction  (x-axis)  instead,  larger  bending  curvature  can  be  achieved. 


Figure  6.  Side  view  and  cross-sectional  view  of  half  sample  of  3-cell  actuator  attached  with  1,  3, 
and  5  segments  of  stiffeners  in  finite  element  simulation  under  applied  nominal  electric  field  of 
22400  V/mm. 

Secondly,  strain  field  from  finite  element  simulation  shown  in  Figure  5  noted  strain  concentration 
around  the  stiffener/elastomer  junctions.  That  is  to  say,  in-plane  expansion  is  not  homogenous 
and  has  localization  related  to  the  location  of  stiffeners.  Each  stiffener  contributes  the  extra 
region  in  addition  to  its  own  width.  Therefore,  multiple  stiffeners  are  able  to  contribute  more 
localization  and  constrain  larger  region  of  in  plane  deformation.  This  gives  us  higher  value  of 
effective  width  that  is  actually  inducing  bending  leads  to  larger  actuation  curvatures. 

Reviewing  analytical  analysis  in  our  previous  work  []  based  on  the  theory  of  bi-material 
thermostats,  axial  strain  is  one  of  the  main  driving  factor  that  urge  bending  and  was  assumed  to 
be  uniform.  Equation  we  used  to  represent  actuation  curvatures  /f  as  a  function  of  input  nominal 
electric  field  was  written  as  equation  (1)  where  ee^  is  dielectric  constant  of  the  elastomer 

multiplies  vacuum  permittivity;  I ,  and  a  are  Young's  modulus,  moment  of  inertia,  and  the 
layer  thickness;  subscripts  1  and  2  represent  layer  stiffener  and  layer  dielectric  elastomer 
respectively,  h  is  the  total  thickness  and  is  a  geometry  factor  that  varies  with  the  amount  of 


laminated  cells  due  to  the  actuator  design  where  in  this  3-cell  case,  A  =  ~  • 
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However,  stiffeners  in  our  design  are  not  fully  covering  the  entire  surface  and,  as  a  result,  only 
approving  an  effective  region  of  constraint.  A  schematic  diagram  in  Figure  7  illustrates  the  idea 
of  partial  axial  deformation  caused  by  single  stiffener  within  an  effective  width.  The  effective 
width  should  be  used  to  modify  our  analytical  analysis.  Furthermore,  the  additional  fitting 
parameter,  0.55,  we  is  found  to  be  caused  by  this  modification.  The  modified  analytical  solution 
is  shown  in  equation  (2)  where  effective  width  Weff  varies  with  different  cases  of  stiffener 
configuration. 
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Figure  7.  Schematic  diagram  of  narrow  strip  stiffener  constrains  only  part  of  the  axial  expansion 
of  EAR.  (a)  In  ideal  case  of  analytical  solution,  uniform  in-plane  expansion  of  EAR  is  expected. 
Total  EAR  width,  \ni  are  affected  by  stiffener,  (b)  In  real  case,  only  partially  constraint  is  caused 
by  narrow  stiffener  reinforcement  and  Weff  is  required  to  replace  W2. 


Fitting  analytical  solution  with  experiential  results,  effective  widths  of  single  and  three  stiffeners 
were  figured  to  be  11mm  and  21mm  respectively.  The  effective  width  can  be  found  through  finite 
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element  simulation  or,  for  this  case, 

'  '  ejj  stijjener  sealing  region 


the  total  width  of  stiffeners,  w  is  the  sealing  region  from  the  actuator  design  which  is 

sealing  region  o  o  o 


3mm  total,  n  \s  the  number  of  stiffeners,  and  is  the  local  constraint  width  caused  by 

each  stiffener,  in  this  case  equal  to  2.5mm.  Results  are  shown  in  Figure  8. 


Figure  8.  Fitting  actuation  curvatures  from  analytical  analysis  \A/ith  experimental  results. 
Actuation  curvatures  shows  deformation  under  applied  nominal  electric  field  of  0-22.4kV/mm  of 
1  (black)  and  3  (red)  segments  of  stiffeners.  Square  marks  and  solid  lines  represent  experiments 
and  finite  element  analysis  respectively. 

Theoretically,  effective  width  can  be  increased  unlimitedly;  however,  after  certain  amount  of 
stiffeners,  the  localization  effective  region  caused  by  each  stiffener  may  become  overlapped.  The 
total  effective  width  is  going  to  be  saturated  with  limited  actuation  presentation.  We  simulated 
the  effect  of  increasing  the  numbers  of  split  stiffeners  from  one  to  five  in  finite  element  model 
as  shown  in  Figure  9Figure  5.  The  increases  of  actuation  curvatures  become  smaller  each  time 
we  split  the  stiffener  into  narrower  segments.  It  saturated  at  four  segments  under  maximum 
applied  electric  field  of  22400  V/mm. 


Figure  9.  Finite  element  results  of  actuation  curvature  versus  applied  nominal  electric  field  of  3- 
cell  planar  actuator  with  different  numbers  of  split  stiffeners  from  one  to  five. 

In  addition,  total  elastic  energy  was  obtained  from  finite  element  analysis  (Figure  10).  The  results 
show  that  the  energy  can  be  considered  as  no  change  with  different  numbers  of  stiffener  we  use. 
The  energy  doesn't  increase  with  higher  actuation  curvatures.  This  tells  us  only  the  efficacy  of 
energy  is  changing.  Actuation  strain  is  more  properly  used  in  the  bending  direction  in  multiple- 
stiffener  case  and  leading  to  higher  bending  curvatures. 
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Figure  10.  Finite  element  results  of  total  elastic  energy  versus  applied  nominal  electric  field  of  3- 
cell  planar  actuator  with  different  numbers  of  split  stiffeners  from  one  to  five. 

Finally,  force-stroke  characteristics  can  also  be  obtained  from  finite  element  simulation.  Take 
applying  electric  field  of  20800  V/mm  as  an  example  which  shown  in  Figure  11,  the  increases  of 
stiffener  periodicity  can  be  used  to  adjust  different  force-stroke  characteristics.  More  stiffeners 
yield  higher  blocking  force  as  well  as  larger  deflection  of  this  bending  device. 
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Figure  11.  Force-stroke  characteristics  of  3-cell  actuator  attached  with  1,  3,  and  5  segments  of 
stiffeners  in  Finite  element  simulation  under  applied  electric  field  of  20800  V/mm  comparing  with 
the  case  of  fully  covered  continuous  stiffener. 


Conclusion 

In  this  paper,  we  investigate  the  role  of  stiffener  in  our  planar  actuator  design.  Stiffeners  bring 


the  ability  out-of-plane  deformation  to  planar  dielectric  elastomer  actuators.  Each  small  segment 
constrains  a  localized  in  plane  expansion  of  elastomer  and  causes  axial  stain  that  leads  the 
actuator  to  bend.  Experimental  and  finite  element  simulation  results  are  provided  to  investigate 
how  changing  stiffener  configuration  while  keeping  the  bending  stiffness  may  influence  actuation 
curvatures.  Moreover,  force-stroke  characteristics  can  be  adjusted  by  varying  stiffener 
periodicity  as  well.  Modified  analytical  analysis  is  also  provided  to  describe  the  effect  of  stiffener 
reinforcement. 


